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ABSTRACT 


A  method  of  designing  airfoils  in  cascade  by  means  of 
conformal  transformations  is  discussed.  Design  of  these  air¬ 
foils  is  regulated  by  five  independent  input  parameters,  with 
solutions  obtained  by  digital  computer.  A  large  number  of 
cascades  generated  by  this  method  are  compared.  To  evaluate 
the  limits  of  performance,  a  parameter  to  indicate  the  ten¬ 
dency  toward  flow  separation  is  introduced,  with  a  limiting 
value  established  and  verified.  Proper  solidity  is  shown  to 
be  of  great  importance  in  achieving  low  values  of  this  separ¬ 
ation  parameter.  The  value  of  proper  solidity  for  a  given 
blade  thickness  is  shown  to  be  relatively  independent  of 
fairing  shape.  To  increase  performance,  reducing  blade  thick¬ 
ness  with  a  corresponding  increase  in  solidity  is  shown  to 
yield  much  greater  improvement  than  changes  in  fairing  shape. 
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1.0  Introduction 

This  study  is  an  extension  of  a  line  of  investigation 
originally  conducted  by  Professor  Theodore  H.  Gawain  of  the 
0.  S.  Naval  Postgraduate  School.  The  original  investigation, 
sponsored  by  the  Aerojet  General  Corporation,  was  specifically 
aimed  at  increasing  the  pressure  rise  per  stage  of  axial  flow 
pumps.  The  intended  application  was  to  reduce  the  size  of 
propellant  pumps  for  large  liquid  fuel  rockets.  The  require¬ 
ments  were  for  maximum  pressure  rise  per  stage,  regardless 
of  efficiency.  This,  in  turn,  demanded  positive  and  acciirate 
control  of  the  pressure  distribution  around  each  blade,  with 
particiilar  emphasis  on  the  adverse  pressure  gradient  experi¬ 
enced  by  the  boundary  layer. 

Flow  of  a  real  fluid  through  an  axial  flow  compressor  or 
pump  involves  some  formidable  problems.  Besides  the  effects 
of  viscosity  and  turbxilence;  tip  losses,  centrifugal  force 
fields,  and  other  effects  complicate  the  problem  considerably. 
As  is  usually  the  case  in  engineering  however,  a  simplifica¬ 
tion  of  the  problem  may  yield  useful  results.  If  the  stream 
surfaces  of  the  flow  are  assumed  to  be  concentric  cylinders, 
then  a  cylindrical  cross  section  of  the  flow  may  be  repre¬ 
sented  by  an  infinite  linear  cascade  with  two-dimensional  flow. 
If  all  velocities  are  considered  relative  to  the  blade  row, 
then  the  solution  for  the  rotor  and  the  stator  have  identical 
mathematical  form.  Following  the  usual  aerod3mamic  practice, 
viscosity  may  be  neglected  outside  the  boundary  layer.  This 
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leaves  only  turbulence  as  a  factor  preventing  the  use  of  poten 
tial  flow  considerations.  In  order  to  use  the  powerful  analy¬ 
tic  methods  of  potential  flow,  turbulence  has  been  neglected. 
This  does  not  mean  however,  that  the  data  presented  is  invalid 
Actually,  the  potential  flow  assumption  is  used  to  obtain  the 
pressure  distribution  only.  The  relation  between  pressure 
distribution  and  separation  was  obtained  from  experimental 
data  from  actual  turbomachines,  and  must  therefore  include  the 
effects  of  turbulence. 

From  the  foregoing  considerations,  it  may  be  seen  that 
the  study  of  cascade  design  has  application  to  the  whole  field 
of  axial  flow  turbomachinery.  The  problem  of  cascade  design 
is  much  more  difficult  than  the  problem  of  isolated  airfoil 
design.  In  the  early  days  of  airfoil  design,  much  experimen¬ 
tal  data  was  obtained  on  a  bewildering  variety  of  shapes. 
However,  there  was  no  reliable  method  of  predicting  the  per¬ 
formance  of  an  airfoil  before  it  was  tested.  It  was  only 
after  the  NACA  (now  NASA)  discovered  how  airfoil  performance 
depended  on  shape,  specifically  on  thickness,  camber  and  fair¬ 
ing  distribution,  that  any  real  progress  was  made  in  the  field 
of  airfoil  design. 

The  present  state  of  the  field  of  cascade  design  is  that 
there  is  a  large  volume  of  experimental  data  available  on  a 
bewildering  variety  of  cascades.  Many  efforts  have  been  made 
to  extrapolate  the  theory  of  isolated  airfoils  so  that  it  will 
apply  to  cascades.  However,  at  present  there  is  no  reliable 
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theory  which  covers  the  field  of  cascade  performance  in  the 
way  that  NASA's  theory  covers  the  field  of  airfoil  performance. 
The  parameters  of  thickness,  camber,  and  fairing  distribution 
are  still  involved,  plus  a  few  more.  The  relations  and  in¬ 
teractions  between  these  parameters  are  so  complex,  that  it 
is  very  difficult  to  consider  any  aspect  of  cascade  design 
separately. 

When  considering  cascades  as  complete  entities,  the  inter 
action  effects  are  all  present,  and  need  not  be  corrected  for. 
However,  experimental  efforts  in  this  direction  have  been  nec¬ 
essarily  limited  in  scope,  due  to  the  effort  and  expense  in¬ 
volved  in  constructing  and  testing  large  numbers  of  cascades. 
Theoretical  methods,  such  as  the  one  outlined  in  this  investi¬ 
gation,  provide  the  means  of  considering  complete  cascades 
over  a  wide  range  of  design  parameters.  This  study  was  under¬ 
taken  from  the  standpoint  of  discovering  areas  of  possible 
improvement  in  compressor  performance. 

Most  of  the  efforts  of  the  NASA  in  the  field  of  compre- 
sor  design  have  been  for  maximtun  efficiency  with  moderate 
output.  In  some  cases,  it  may  be  required  to  attain  maximum 
output,  with  perhaps  moderate  efficiency.  The  gains  to  be 
expected  from  these  considerations  are  not  of  large  order  how¬ 
ever,  as  a  cascade  develops  its  best  efficiency  not  far  from 
the  point  of  stalling  the  blades. 

The  basic  problem,  then,  is  composed  of  two  parts.  First 
it  is  necessary  to  devise  a  method  for  generating  blade  shapes 
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vd.th  the  corresponding  pressure  distributions.  For  this  pur¬ 
pose,  the  powerful  methods  of  complex  variables  and  conformal 
transformations  are  utilized.  The  method  developed  by  Pro¬ 
fessor  Gawain  generates  infinite  cascades  of  airfoil  shapes, 
based  on  five  arbitrary  parameters.  Generally,  the  transfor¬ 
mations  are  similar  to  the  Joukowski  transformations  with 
additional  requirements. 

The  second  part  of  the  problem  is  to  establish  a  correla¬ 
tion  between  the  pressure  distribution  around  each  blade  and 
the  occurence  of  separation.  This  problem  must,  of  course, 
depend  heavily  on  experimental  data.  Fortunately,  work  has 
been  done  in  this  area  and  some  practical  working  limits  have 
been  established.  See  Reference  1.  Based  on  this  data,  the 
simple  parameter  (  been  chosen  as  an  indication 

of  the  tendency  toward  separation. 

The  author  wishes  to  express  his  sincere  appreciation  to 
Professor  Theodore  H.  Gawain  of  the  U»  S.  Naval  Postgraduate 
School  for  his  assistance,  guidance,  and  encouragement  during 
the  course  of  this  investigation. 
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2,0  The  Method  of  Cascade  Design 


2.1  Introduction 

Methods  of  designing  airfoils  by  purely  mathematical  means 
have  been  known  since  the  early  days  of  aviation,  Joukowski 
was  the  most  notable  pioneer  in  this  field.  He,  and  others, 
noted  the  possibility  of  designing  cascades  of  airfoils  by 
similar  methods.  Since  the  early  applications  of  cascade  de¬ 
sign  were  to  multi-winged  airplanes,  not  much  work  was  done  in 
this  field.  As  pointed  out  earlier,  the  application  to  turbo¬ 
machinery  makes  this  field  of  considerable  interest  today. 

In  considering  the  flow  through  cascades,  the  shape  and 
position  of  any  Individual  blade  will  affect  the  flow  around 
that  blade.  In  addition,  the  shape  and  position  of  all  the 
other  blades  will  also  affect  the  flow  around  that  same  blade 
at  the  same  time.  Therefore,  more  parameters  must  be  used  in 
cascade  design  than  are  used  in  isolated  airfoil  design.  In 
addition  to  blade  shape,  the  relation  of  the  blades  to  each 
other,  and  their  relation  to  the  flow  must  all  be  described. 

The  method  developed  by  Professor  Gawain  uses  five  independent 
parameters  to  describe  an  individual  cascade. 

The  solution  for  blade  shape  and  velocity  distribution 
corresponding  to  chosen  values  of  the  five  parameters  is  found 
by  using  relations  of  potential  flow  and  conformal  mapping  in 
three  separate  planes.  These  have  been  termed  the  near  circle 
plane,  the  circle  plane,  and  the  cascade  plane  respectively. 
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In  reading  through  the  follo'fcdng  sections,  which  describe 
the  general  method  of  generating  cascades,  reference  may  be 
made  to  Figure  5t  which  illustrates  a  typical  contour  in  each 
of  the  planes  referred  to  above.  A  summary  of  the  mathematical 
development  of  the  method  is  given  in  Appendix  A.  A  more  de¬ 
tailed  analysis  may  be  foxmd  in  Reference  2. 

The  solution  of  these  equations,  many  of  which  are  trans¬ 
cendental,  is  an  impossible  task  by  hand  calculation.  The 
entire  problem  has  been  arranged  in  FORTRAN  language  for  solu¬ 
tion  by  a  digital  computer.  The  GDC  1604  computer  at  the  U.  S. 
Naval  Postgraduate  School  was  used  for  this  investigation. 

The  basic  computer  program,  along  with  later  modifications  is 
included  in  Appendix  B. 

The  following  sections  describe  the  general  method  of 
generating  cascades  from  five  input  parameters.  Each  parameter 
is  described  in  detail  as  it  enters  into  the  calculations. 

The  parameters  are; 

=  Mean  flow  direction 
do  “  Turning  parameter 
0^  -  Equivalent  flat  plate  solidity 

S'  “  Thickness  parameter 
/<  =  Shape  or  symmetry  parameter 

2.2  The  Flow  Parameters 

Flow  through  a  cascade  is  usually  described  in  terms  of 
two  velocities.  The  first  velocity  is  that  far  enough  up- 
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stream  of  the  cascade  so  that  local  effects  due  to  the  cascade 
are  not  evident.  This  is  called  the  inlet  velocity,  and  the 
angle  it  makes  with  the  reference  axis  is  called  /3j^.  Similarly, 
the  velocity  far  enough  downstream  so  that  local  effects  have 
damped  out  is  called  the  outlet  velocity,  and  its  angle  with 
the  reference  axis  is  called  ^3.  The  magnitudes  of  these  two 
velocities  are  not  independent,  since  the  laws  of  continuity 
require  that  the  component  perpendicular  to  the  cascade  axis 
be  equal  for  the  two  cases.  This  concept  allows  velocities  to 
be  conveniently  nondimensionalized,  so  that  only  angles  are 
important . 

Figure  1  is  an  illustration  of  two  typical  cascades.  An 
appropriate  velocity  vector  diagram  is  shown  for  each.  Note 
that  only  two  parameters  are  necessary  to  describe  the  complete 
diagram.  Instead  of  the  angles  (3i  and  ^33  it  was  more  conven¬ 
ient  to  use  two  different  parameters.  Consider  the  following 
expressions: 

Tan  /$  “  i(Tan  +  Tan  /32) 

Bq  =  (Tan  pi  -Tan  (32) 

The  first  expression  defines  the  angle  p,  or  mean  flow  direc¬ 
tion,  which  is  the  first  input  parameter.  The  second  expres¬ 
sion  defines  the  turning  parameter  Bq,  which  is  the  second 
input  parameter.  The  relationship  between  all  these  quantities 
is  shown  in  the  diagram  on  Figure  1. 

The  starting  point  of  the  calculations  is  the  near  circle 
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plane.  In  this  plane ,  the  complex  potential  for  a  series  of 
doublets  spaced  at  intervals  of  2  JT  along  the  y  axis  is  com¬ 
bined  with  the  potential  for  uniform  flow  at  the  angle  with 
the  X  axis.  The  velocity  of  this  flow  is  set  at  unity ,  which 
then  becomes  the  reference  for  all  the  nondimensional  veloci¬ 
ties  derived  later.  If  only  one  doublet  were  used,  the  resul¬ 
ting  potential  would  describe  the  flow  without  circulation 
around  a  single  circular  cylinder.  In  this  case,  with  an  in¬ 
finite  series  of  doublets  along  the  y  axis,  the  resulting 
potential  describes  the  flow  without  circulation  around  an 
infinite  series  of  nearly  circular  bodies  spaced  along  the  y 
axis.  These  bodies  are  not  actually  circles  because  the  pre¬ 
sence  of  the  infinite  stack  of  doublets  distorts  the  flow 
around  any  individual  doublet,  such  that  the  flow  contour  is 
a  slightly  flattened  circle.  The  amount  of  distortion  is  re¬ 
lated  to  the  strength  of  the  doublets  and  the  spacing  between 
them.  In  this  case,  the  spacing  between  doublets  is  fixed  at 
2  77*  j  so  doublet  strength  is  vrtiat  determines  the  relation 
between  spacing  and  size  of  the  bodies  in  the  near  circle  plane. 

2.3  The  Solidity  Parameter 

The  strength  of  the  doublets  is  fixed  by  the  input  para¬ 
meter  <y  ,  or  the  equivalent  flat  plate  solidity.  This  para¬ 
meter  is  closely  related  to,  but  not  identical  with, the  ordinary 
solidity  c/S,  or  the  chord-pitch  ratio  of  the  blade  row.  The 
distinction  can  best  be  illustrated  by  an  example.  Generally, 


any  row  of  arbitrary  blades  of  solidity  c/S  may  be  transformed 
conformally  into  a  related  row  of  flat  plate  vanes.  The  chord- 
pitch  ratio  of  the  latter  is  then  denoted  by  G'  .  It  follows 
that  if  two  arbitrary  cascades  have  the  same  value  of  , 
then  either  may  be  conformally  transformed  into  the  other.  But 
if  the  value  of  CT"  is  different  for  each,  then  no  such  trans¬ 
formation  is  possible.  then,  may  be  regarded  as  a  kind 

of  generalized  solidity,  of  perhaps  deeper  aerod3mamic  signi¬ 
ficance  than  simple  geometric  solidity. 

The  relation  between  values  of  CT" and  geometric  solidity 
is  directly  influenced  by  blade  thickness.  Very  thin  blades 
have  geometric  solidities  almost  exactly  equal  to  CT* .  With 
thicker  blades,  the  geometric  solidity  is  always  less.  The 
reason  for  this  may  be  illustrated  by  examples  of  two  extremes. 
If  a  contour  in  the  near  circle  plane  is  transformed  confor¬ 
mally  into  a  flat  plate,  the  chord  length  of  the  flat  plate 
will  be  nearly  twice  the  diameter  of  the  original  contotir. 
However,  if  no  transformation  is  made  at  all,  an  "airfoil"  of 
100^  thickness  results,  whose  chord  length  is  equal  to  the 
diameter.  Both  cascades  would  have  the  same  value  of  G^,  but 
the  geometric  solidity  of  the  latter  case  would  be  half  that 
of  the  former.  Between  these  two  extremes,  finite  thickness 
airfoils  in  cascade  have  solidities  less  than  . 

2.4  The  Turning  Parameter 

The  turning  parameter,  Bq,  is  a  measure  of  the  change  in 
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direction  of  the  flow  passing  thi’otigh  the  blade  row.  It  is 
actually  a  nondimensional  form  of  the  downwash  velocity  induced 
by  the  blades.  There  is  a  very  simple  relation  between  Bq  and 
the  circulation  around  each  blade.  In  fact,  for  a  constant 
mean  flow  angle,  this  relation  is  a  simple  direct  proportion. 
This  convenient  state  of  affairs  is  the  reason  why  the  two 
flow  parameters  were  expressed  in  their  present  form. 

In  isolated  airfoil  theory,  the  method  of  introducing  lift 
or  circulation  on  an  airfoil,  is  to  add  to  the  complex  poten¬ 
tial,  a  term  representing  a  vortex  of  suitable  strength.  How¬ 
ever,  mere  addition  of  a  vortex  will  produce  a  distortion  in 
the  contours  in  the  near  circle  plane,  since  these  contoxirs 
are  not  exactly  circles.  In  addition,  an  infinite  stack  of 
vortices  is  necessary  to  produce  circulation  aroiind  each  of  the 
blades  in  the  cascade.  To  avoid  these  resulting  undesirable 
pertxarbations  in  what  will  ultimately  become  the  airfoil  shape, 
and  indirect  method  of  adding  circulation  is  used. 

In  this  method,  a  transformation  is  made  from  the  near 
circle  plane  to  the  circle  plane.  An  equation  is  developed 
which  transforms  the  flow  in  the  near  circle  plane  into  the 
flow  around  a  perfect  circle,  centered  on  the  origin.  The  cir¬ 
culation  term  may  then  be  introduced  without  fear  of  distor¬ 
tion.  The  transformation  is  then  made  back  to  the  near  circle 
plane.  Since  the  same  equation  is  used  to  transform  both  ways, 
the  original  flow  contour  remains  unchanged;  however,  the 
velocities  at  each  point  are  now  different,  due  to  the  circu- 
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lation 


2,5  The  Blade  Shape  Parameters 

The  final  transformation  is  from  the  near  circle  plane 
to  the  cascade  plane.  This  is  done  by  the  Joukowski  method 
of  shifting  the  axes  slightly  in  such  a  way  that  the  flow  con¬ 
tour  transforms  into  an  airfoil  shape.  Two  arbitrary  para¬ 
meters  were  chosen  to  control  the  blade  shape.  One  of  these 
is  called  the  thickness  parameter,  and  is  given  the  symbol  (S  . 
This  quantity  is  generally  related  to  the  thickness,  and  also 
to  nose  curvature,  but  not  in  any  fixed  relationship.  Large 
values  of  £,  produce  thick  airfoils,  and  small  values  produce 
thin  airfoils.  Figure  2  is  an  illustration  of  several  airfoils, 
showing  the  effect  of  .  The  limiting  value  of  0  for  &  does 
not  however  produce  an  airfoil  of  zero  thickness.  Due  to  the 
slight  distortion  of  the  contours  in  the  near  circle  plane, 

=  0  may  produce  the  mathematically  significant,  but  physi¬ 
cally  useless  case  where  the  "lower"  surface  of  the  airfoil 
crosses  throixgh  the  "upper"  surface,  then,  can  be  small  but 
not  zero,  and  the  airfoils  produced  can  be  thin  but  not  of  zero 
thickness. 

The  other  parameter  is  the  shape  or  symmetry  parameter  K, 
The  effect  of  variations  in  K  may  be  seen  in  Figure  3®  Gener¬ 
ally,  K“0  produces  a  blunt  leading  edge,  and  cusped  trailing 
edge  similar  to  a  Jotikowski  airfoil,  K=1  will  produce  a  sym¬ 
metrical  shape  which  is  generally  a  cambered  ellipse.  As  K 
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increases  from  zero  to  one,  the  point  of  maximum  thickness 
moves  from  the  quarter  chord  to  the  half  chord  point,  and  the 
trailing  edge  roxinds  offo  Values  of  K  larger  than  one  will 
produce  airfoils  facing  the  wrong  way  (Trailing  edge  into  the 
flow)  and  hence  must  not  be  usedo 

2.6  Stagnation  Points 

During  the  transformations  made  in  this  analysis,  the 
stagnation  points  serve  as  reference  points  to  locate  the  re¬ 
lation  of  the  contours  in  the  various  planes.  For  this  reason 
there  is  always  a  stagnation  point  at  the  leading  and  the 
trailing  edge.  Having  a  stagnation  point  at  the  leading  edge 
is  sometimes  called  the  Theodorsen  condition.  This  means  that 
in  all  the  cascades  developed  by  this  method,  the  airfoils  are 
operating  at  the  ideal  angle  of  attack.  Variations  in  lift 
then,  are  achieved  by  variations  in  the  camber  of  the  blades. 
Figure  1  graphically  illustrates  this  fact.  Two  typical 
cascades  are  drawn,  with  the  principal  variation  being  in  turn¬ 
ing  parameter,  The  values  of  ^  were  adjusted  to  give 

the  same  blade  thickness.  The  difference  in  camber  may  be 
clearly  seen. 

The  condition  of  a  stagnation  point  at  the  trailing  edge 
is  called  the  Kutta  condition.  This  fact  has  particular  sig¬ 
nificance  when  considering  the  blade  shapes  where  K=l.  As 
pointed  out  previously,  these  shapes  are  generally  cambered 
ellipses.  It  is  a  well  known  fact  that  flow  around  an  ellipse 
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of  a  given  thickness  ratio  produces  a  lower  velocity  peak  than 
the  flow  around  any  other  shape  of  the  same  thickness  ratio. 
This  class  of  blades  then  should  exhibit  the  lowest  loading 
attainable  for  any  possible  fairing  shape.  Of  course,  any 
attempt  to  duplicate  this  fairing  shape  in  an  actual  cascade 
is  doomed  to  failure  because  of  the  rounded  trailing  edge. 
However,  this  case  still  has  significance  in  two  respects. 
First,  the  performance  of  a  cambered  ellipse,  with  the  Kutta 
condition  mathematically  imposed,  represents  a  useful  limiting 
case  of  performance  due  to  fairing  shape.  Secondly,  a  blade 
similar  to  this,  with  modifications  to  produce  a  sharp  trail¬ 
ing  edge,  would  have  very  nearly  the  same  perfonnance.  The 
critical  area  of  flow  near  the  leading  edge  would  not  be  af¬ 
fected  much  by  modifications  at  the  rear  of  the  airfoil. 

Figure  4  is  a  plot  of  the  pressure  distribution  around 
the  airfoils  of  Figure  3*  The  above  premise  is  verified,  as 
the  negative  pressure  peak  can  be  seen  to  decrease  as  the 
value  of  K  is  increased.  Note  that  the  pressture  distribution 
aft  of  the  peak  is  essentially  linear  to  the  trailing  edge. 
Reference  3,  p.  354  states  that  this  type  of  pressure  distri¬ 
bution  is  typical  of  decelerating  cascades.  Note  also  that 
the  slope  of  this  line  becomes  less  as  K  is  increased,  indi¬ 
cating  a  less  adverse  pressure  gradient  in  the  boundary  layer. 
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3.0  Prediction  of  Separation 


3.1  Pressure  Gradient 

The  phenomenon  which  ultimately  limits  the  attainable 
pressure  rise  of  a  cascade  is  separation.  Flow  separation 
occurs  vrtien  flow  in  a  boundary  layer  experiences  an  adverse 
(positive)  pressure  gradient.  This  adverse  gradient  tends  to 
decrease  the  velocities  in  the  boundary  layer.  Ultimately, 
the  point  is  reached  where  the  velocities  close  to  the  surface 
have  been  reduced  to  zero.  At  this  point,  flow  separation 
occurs.  Beyond  this  point,  some  velocities  become  negative, 
and  a  region  of  reverse  flow  exists.  The  resultant  effects 
cause  high  drag  and  loss  of  circulation  around  each  blade 
element.  Many  experiments  have  verified  that  any  adverse 
pressure  gradient,  however  small,  will  cause  separation  if  it 
acts  long  enough.  The  concept  of  a  gradient  and  the  distance 
over  which  it  acts,  leads  one  naturally  to  suspect  pressure 
d4.1'ierence  as  a  significant  parameter  in  predicting  separation. 
The  gradient,  integrated  over  the  distance  is  mathematically 
equal  to  the  difference  in  pressure.  It  was  previously  shown, 
furthermore,  that  the  pressure  distribution  over  typical  air¬ 
foils  in  this  system  is  approximately  linear  aft  of  the  nega¬ 
tive  pressure  peak, 

3.2  Velocity  Considerations 

From  another  point  of  view,  one  might  examine  the  velo¬ 
cities  of  flow  through  a  cascade.  From  the  inlet  velocity 
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described  previously,  velocity  is  brought  to  zero  at  the 
front  stagnation  point,  accelerated  to  a  local  maximum  some¬ 
where  on  the  surface  of  the  blades,  reduced  toward  zero  at 
the  trailing  edge,  and  finally  settles  down  to  the  outlet 
velocity  far  behind  the  cascade.  Besides  the  inlet  and  outlet 
velocities,  only  one  other  velocity  is  really  significant,  and 
that  is  the  local  maximum,  .  From  these  three  quantities, 

only  two  independent  dimensionless  parameters  may  be  construc¬ 
ted.  If  the  inlet  velocity  and  were  used,  this  parameter 

would  be  related  to  the  flow  aroxmd  the  forward  part  of  the 
airfoil.  Since  this  flow  involves  a  negative  pressure  grad¬ 
ient,  it  is  probably  not  very  significant  with  respect  to 
separation.  Similarly,  the  quantities  V___  and  the  outlet 
velocity,  called  V2,  are  related  to  the  flow  around  the  rear 
part  of  the  airfoil,  which  involves  a  positive  pressure  grad¬ 
ient.  These  two  parameters  then,  combined  in  some  dimension¬ 
less  manner,  should  act  as  an  indicator  for  the  tendency 
toward  separation. 

3.3  Diffusion  Factors 

From  its  many  experimental  studies  of  flow  through  cas¬ 
cades,  the  NASA  has  arrived  at  essentially  the  same  conclusion. 
They  have  used  several  dimensionless  combinations  of  and 

V2  to  define  blade  loading  limits.  One  of  these  is  the 
’’Local  Diffusion  Factor”  which  is: 
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(7  -  V  ) 

D,  “  max  2 
V 

max 

This  factor  is  described  in  Reference  1  and  supporting  data 
indicates  that  separation  is  likely  to  occur  for  values  of 
greater  than  0«5« 

Since  the  velocity  distribution  around  the  blades  must 
be  known  to  find  values  of  ,  it  is  necessary  to  use  experi¬ 
mental  data,  or  an  analytical  treatment  such  as  the  one  des¬ 
cribed  in  this  investigation.  To  get  aroxind  this  difficulty, 
the  NASA  also  defines  an  approximation,  D,  which  may  be  ob¬ 
tained  from  more  easily  determined  quantities.  D,  is  called 
the  ’’Diffusion  Factor"  and  is  given  by: 

D  **  1-  Cos  +  Cos  (Tan/5l  -  Tany52) 

Cos  (^2  Vs 

It  is  seen  that  this  quantity  can  be  determined  from  the  inlet 
and  outlet  velocities,  and  the  solidity  of  the  blade  row.  The 
supporting  data  indicates  that  separation  is  likely  to  occxir 
for  values  of  D  greater  than  0.6. 

Since  the  methods  used  in  this  investigation  furnish  all 
of  the  flow  quantities  described  above,  as  well  as  the  com¬ 
plete  velocity  distribution  around  each  blade,  the  opportunity 
to  compare  the  two  NASA  parameters  is  presented.  As  will  be 
described  later,  the  correlation  between  these  two  parameters 
is  not  too  good  for  the  solidities  considered  here. 
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3«4  The  Boundary  Layer  Loading  Parameter 

From  the  above  development,  it  may  be  seen  that  the  para¬ 
meter  to  indicate  separation  should  involve  the  pressxire  ratio 
across  the  rear  part  of  the  airfoil,  and  be  some  dimensionless 
combination  of  the  quantities  V2  and  The  NASA  has  in¬ 

dicated  that  most  compressors  operate  in  the  range  where 
Reynolds  number  effects  are  insignificant.  Therefore,  it  is 
felt  that  no  particular  inaccuracies  result  from  omitting 
consideration  of  Reynolds  number. 

In  order  to  avoid  using  one  parameter  for  pressure  dis¬ 
tribution,  and  another  for  the  tendency  toward  separation, 
the  quantity  S  was  selected  vrtiere; 

This  quantity  is  referred  to  as  the  boundary  layer  loading 
parameter,  and  is  quite  versatile  in  application.  It  is 
directly  related  to  the  pressure  coefficient,  and  in  fact  only 
differs  from  that  quantity  by  unity,  S  also  has  a  constant 
functional  relationship  with  the  NASA  parameter  D^o  The  value 
S  =  4*25  corresponds  to  the  value  of  of  0,5,  and  has  been 
selected  as  the  limit  loading  in  the  boundary  layer  for  blades 
in  this  investigation. 
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4.0  Results  and  Discussion 

Theory  and  experiment  both  point  toward  increasing 
solidity  as  the  method  of  obtaining  increased  pressure  rise 
per  stage.  As  solidity  is  increased,  the  lift  coefficient  on 
each  blade  is  decreased,  and  hence  the  loading  on  the  boundary 
layer  is  decreased  also.  However,  as  the  solidity  is  increased, 
there  is  an  adverse  effect  due  to  the  blockage  caused  by  the 
finite  thickness  of  the  blades  plus  the  boundary  layer o 
Finally,  the  point  is  reached  where  the  beneficial  effects  of 
decreased  blade  loading  are  just  counterbalanced  by  the  ad¬ 
verse  effects  of  blockage.  The  quantitive  determination  of 
this  point  was  one  of  the  major  objectives  of  this  investiga¬ 
tion. 


4.1  Organization  and  Scope 

Since  actual  cascades  are  designed  on  the  basis  of  blade 
thickness  and  geometric  solidity,  it  was  felt  that  the  results 
of  this  investigation  ought  to  be  organized  on  this  basis, 
rather  than  on  the  basis  of  the  mathematical  parameters  and 
.  The  value  of  blade  thickness  is  not  directly  predictable 
from  the  input  parameter  C  .  Therefore ,  the  basic  computer 
program  was  modified  to  use  an  iterative  scheme  to  achieve  the 
desired  values  of  blade  thickness.  Details  are  given  in 
Appendix  B,  which  describes  the  basic  computer  program  and  the 
modifications.  Seven  blade  thicknesses  were  selected  to  cover 
the  field  of  possible  blades.  These  ranged  from  the  limiting 
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case  of  extremely  thin  blades  to  ithickness  blades  in  four 
percent  increments* 

It  was  decided  to  conduct  this  investigation  using  only 
one  value  of  average  flow  angle  /S*  This  value  was  chosen  to 
be  45 which  was  a  somewhat  arbitrary  decision,  but  not  with¬ 
out  some  sound  theoretical  basis.  It  is  stated  in  Reference 
4  that  maximum  efficiency  may  be  obtained  when  the  average 
flow  angle  is  near  45®,  and  many  compressors  are  built  accord¬ 
ingly. 

In  spite  of  reducing  the  parameters  to  four  by  the  above 
decision,  over  1470  separate  cascades  were  generated  and  ana¬ 
lyzed  for  this  investigation.  Additional  parameters  could 
have  been  introduced  to  give  more  control  over  the  blade  shape, 
but  they  would  have  greatly  magnified  the  problem.  As  the 
data  will  show,  additional  blade  shape  parameters  would  not 
have  affected  the  significance  of  the  results  very  much. 

4.2  Series  1  Plots 

The  first  plots  of  data  from  the  computer  run  were  made 
of  boundary  layer  loading  vs.  solidity,  for  airfoils  of  con¬ 
stant  blade  thickness.  Five  values  of  the  shape  parameter  K 
were  used,  and  data  was  obtained  at  six  values  of  turning 
parameter,  This  meant  that  there  were  30  graphs  in  series 

one.  Figiu’es  6  through  9  are  representative  graphs  from  this 
series. 

It  was  noted  with  satisfaction  that  the  curves  showed  a 
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pronounced  minimum  loading  as  soliaity  was  increased.  The  one 
exception  was  the  case  of  the  extremely  thin  airfoils,  for 
which  the  curve  appeared  to  be  asymptotic  to  some  low  value  of 
loading*  This  is  not  unreasonable,  since  the  blockage  effect 
would  be  very  small  for  this  case.  The  case  of  thin  airfoils 
then,  represent5a  sort  of  limiting  value  for  cascades.  The 
point  at  which  minimum  loading  occuired  for  a  given  thickness 
airfoil  was  called  the  optimtim  solidity. 

It  may  be  noted  that  while  the  value  for  loading  at  an 
optimum  point  could  be  determined  with  accijracy,  the  corres¬ 
ponding  solidity  for  that  point  was  subject  to  some  uncertainty. 
Therefore,  another  point  was  introduced  where  the  solidity 
could  be  determined  with  greater  certainty.  This  was  called 
the  design  point,  and  is  defined  as  the  point  toward  lower 
solidity  where  the  loading  equals  1.10  times  the  minimum  value. 
There  is  further  justification  in  this  step  from  the  fact  that 
the  basic  analysis  neglects  boundary  layer.  The  blockage 
effect  of  a  blade  plus  boundary  layer  must  be  greater  than  the 
effect  due  to  the  blade  alone.  Therefore,  the  true  optimum 
solidity  must  be  somewhat  less  than  the  apparent  optimiim  based 
on  non-viscous  flow.  The  exact  relation  cannot  be  found  with¬ 
out  a  detailed  analysis  of  the  boundary  layer.  This  simple 
method  sidesteps  this  formidable  problem,  and  is  at  least 
qualitatively  in  the  right  direction. 

For  the  case  of  higher  values  of  turning  parameter  and 
thick  blades,  there  was  some  scatter  in  the  data.  Since  the 
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data  is  all  analytical,  and  shoxild  therefore  be  reasonably 
regular,  some  effort  was  made  to  find  the  source  of  the  scat¬ 
ter.  The  relationship  between  the  parameters  <J'  and  <^/S  was 
suspected  first.  Comparative  plots  of  these  two  parameters 
however,  showed  that  they  varied  in  a  regular  fashion  with 
blade  thickness.  It  was  then  decided  that  the  scatter  was  due 
to  the  fact  that  velocities  were  meastired  only  at  40  points 
along  the  airfoil.  If  the  true  maximum  velocity  occixrred  be¬ 
tween  these  finite  points,  then  the  value  for  boundary  layer 
loading  parameter  would  be  in  error.  Remembering  the  fact 
that  this  error  could  be  only  on  the  low  side  was  an  aid  in 
fairing  the  ciurves  through  the  scattered  points. 

It  was  noted  that  the  data  became  more  regular  as  the 
value  of  K  approached  one.  For  these  airfoils,  either  the 
flow  was  more  regular  about  the  nose,  with  no  sharp  peaks  in 
velocity,  or  the  point  of  maximum  velocity  happened  to  coin¬ 
cide  more  often  with  one  of  the  finite  points.  Examination 
of  the  pressure  plots  of  Figure  4  points  toward  the  former 
as  the  cause.  With  high  values  of  K,  the  pressure  changes  are 
seen  to  be  more  gradual,  and  of  lower  magnitude.  Because  of 
this,  the  actual  maximum  velocity  was  always  close  to  the 
velocity  recorded  at  one  of  the  40  points. 

Figure  10  is  an  illustration  of  the  effect  of  K  in  reduc¬ 
ing  the  boundary  layer  loading  for  a  given  value  of  turning 
parameter.  With  all  other  factors  equal,  the  loading  steadily 
decreased  as  the  value  of  K  increased  toward  unity.  Note  that 
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the  effect  of  K  was  much  greatex*  lox"  the  thick  blades.  For 
the  case  of  extremely  thin  blades,  K  had  no  effect  at  all, 
which  was  not  an  unexpected  result. 

4«3  Series  2  Plots 

Since  the  original  series  involved  30  graphs,  some  effort 
was  made  to  summarize  this  data.  The  result  of  this  effort 
was  Figures  11  through  15,  which  collect  all  the  data  for  a 
particular  value  of  K  on  a  single  page.  Boundary  layer  load¬ 
ing  parameter  and  solidity  are  plotted  along  the  coordinates 
as  before,  but  now  the  solidity  is  design  solidity.  In  effect, 
all  the  design  points  for  a  particular  value  of  K  have  been 
plotted.  Points  of  constant  values  of  turning  parameter  have 
been  connected,  as  have  points  of  constant  airfoil  thickness. 
The  result  is  that  four  interdependent  parameters  are  shown 
for  each  point  on  the  chart.  Any  chosen  point  then,  repre¬ 
sents  a  particular  value  of  blade  thickness  and  turning  para¬ 
meter,  with  the  corresponding  botindary  layer  loading  and 
design  solidity  for  the  cascade. 

Figures  16  through  20  are  plots  made  in  exactly  the  same 
manner,  except  that  optimum  points  were  used  instead  of  design 
points.  It  was  necessary  to  change  the  scale  for  solidity 
along  the  horizontal  axis,  since  optimum  points  represent 
higher  values  of  solidity  in  every  case.  In  addition  some 
difficulty  was  experienced  with  irregularity  in  the  data. 

This  was  due  at  least  in  part  to  the  uncertainty  in  determin- 
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ing  the  value  of  solidity  at  the  optimum  point o  This  optimtim 
point  data  was  included  as  a  liraitiing  case,  to  illustrate  the 
maximum  solidities  to  be  used  when  the  boimdary  layer  is  neg¬ 
ligibly  thin. 

The  choice  of  solidity  as  a  coordinate  is  slightly  mis¬ 
leading,  It  might  be  inferred  from  this  choice  that  solidity 
should  be  used  as  a  design  input.  However,  solidity  was 
plotted  in  this  manner  purely  as  a  convenience  in  transferring 
data  from  the  series  one  graphs.  In  design  work,  the  two 
"inputs”  could  be  boundary  layer  loading,  and  turning  para¬ 
meter.  The  "output"  information  would  then  be  blade  thickness 
and  solidity.  An  attempt  was  made  to  plot  boundary  layer 
loading  parameter  vs,  turning  parameter,  but  the  resulting 
blade  thickness  and  solidity  lines  were  so  nearly  parallel 
that  the  information  was  difficult  to  read, 

4.4  Comparison  of  NASA  Separation  Parameters 

The  data  on  the  series  2  plots  contains  all  of  the  infor¬ 
mation  necessary  to  compute  both  NASA  separation  parameters. 
This  affords  an  opportunity  to  compare  these  two  parameters 
directly.  The  boundary  layer  loading  parameter  S,  and  the 
NASA  parameter  have  a  constant  relationship  such  that: 

Thus  a  equal  to  0,5  corresponds  to  a  value  of  S  of  4o25o 
This  value  is  entered  as  a  heavy  line  on  all  series  2  plots. 


23 


The  approximation  D,  is  relacea  to  the  parameters  B, 
and  c/S  by  the  following  relationship; 


Thus,  since has  been  fixed  at  one  value,  each  value  of  Bq 
and  c/S  will  result  in  a  particular  value  for  Do  This  expres¬ 
sion  for  D  shows  that  it  cannot  be  a  good  approximation  over 
a  wide  range  of  solidities,  since  it  indicates  that  loading 
will  always  decrease  whenever  solidity  is  increasedo  It  has 
already  been  demonstrated  that  as  solidity  increases  beyond 
the  optimum  point,  loading  increases  also  due  to  the  blockage 
effect. 

Lines  corresponding  to  a  value  of  D=o5  and  D*=o6  have  been 
plotted  on  Figures  11  through  20  as  dashed  lines.  A  note  of 
caution  is  necessary  in  using  the  data  on  these  figures.  It 
must  be  remembered  that  there  is  really  only  one  solidity 
shown  on  each  chart,  and  that  is  the  design  (or  optimum)  soli¬ 
dity.  The  fact  that  solidity  is  shown  as  a  variable  is  due 
solely  to  the  fact  that  each  combination  of  blade  thickness 
and  turning  parameter  demands  a  particular  value  of  solidity 
at  the  design  (or  optimum)  point. 

If  the  parameter  D  is  a  true  approximation  for  then 
the  lines  corresponding  to  D-.5  and  Di^.5  should  be  nearly 
the  same.  Examination  of  the  figures  shows  that  there  is  a 
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fair  match  for  the  case  of  optimum  solidity.  Generally,  D 
gives  a  value  of  loading  which  is  higher  than  D^,  but  the 
error  is  only  around  10^,  which  is  not  unreasonable.  For  the 
case  of  design  solidity,  there  is  a  greater  discrepancy.  For 
very  thick  blades,  D  still  indicates  a  value  of  loading  higher 
than  D^,  but  this  relation  reverses  so  that  for  thin  blades 
indicates  the  higher  loading. 

The  data  in  Reference  1  states  that  separation  is  likely 
to  occur  for  values  of  greater  than  ,5  or  for  values  of  D 
greater  than  .6,  This  in  itself  is  an  indication  of  a  mis¬ 
match  between  D  and  D^,  If  the  occurrence  of  separation  is 
taken  as  the  matching  criterion,  then  the  line  corresponding 
to  D=.6  should  match  the  D^=.5  line.  For  the  case  of  optimum 
solidity  it  may  be  seen  that  the  matching  is  very  poor,  where 
D  indicates  a  much  higher  loading  than  At  design  solidity 

the  matching  is  better,  but  D  still  indicates  values  of  load¬ 
ing  which  are  around  IS^  higher  than  Extrapolating  this 

trend  leads  one  to  suspect  that,  for  some  value  of  solidity 
less  than  design  solidity,  the  two  parameters  D=c6  and  D3^=,5 
would  probably  match  quite  closely.  From  this,  it  is  not  un¬ 
reasonable  to  expect  that  the  NASA  obtained  its  experimental 
data  from  cascades  whose  solidity  was  less  than  the  design  or 
optimum  as  defined  in  this  report,  A  few  checks  of  source 
data  for  Reference  1  have  verified  this  premise.  Solidities 
were  generally  1.0  for  cascades  with  blade  thicknesses  of 
eight  and  ten  percent,  A  check  of  the  series  2  plots  will 
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reveal  that  this  is  less  than  the  assign  or  optimxira  point  for 
any  fairing  shape. 

4*5  Comparison  of  Blade  Thickness  and  Fairing  Shape  Effects 

Figures  11  through  20  clearly  show  the  gains  that  can  be 
made  in  turning  parameter,  at  no  increase  in  loading,  by 
using  thinner  blades  at  the  corresponding  higher  solidities. 
However,  it  was  desired  to  compare  the  effects  of  changes  in 
fairing  shape  and  blade  thickness  on  the  same  chart.  Figure 
21  was  made  using  the  design  point  data,  and  Figure  22  was 
made  with  data  from  the  optimum  points.  The  value  of  S  =  4*25 
was  selected  as  the  best  estimate  of  maximum  permissible  blade 
loading.  The  values  plotted  were  turning  parameter  vs.  solid¬ 
ity  at  this  constant  value  of  blade  loading.  This  resiilted  in 
lines  which  represented  varying  blade  thickness  at  a  constant 
value  of  K.  Lines  of  constant  blade  thickness  were  then  drawn 
in.  It  is  interesting  to  note  that  these  lines  are  nearly 
vertical.  This  means  that  the  value  of  design  or  optimum 
solidity  for  a  particular  blade  thickness  is  relatively  inde¬ 
pendent  of  blade  shape.  This  result  was  not  directly  predic¬ 
table  from  theory,  and  may  prove  useful  in  planning  experi¬ 
mental  investigations. 

As  blade  thickness  is  decreased,  it  is  clearly  seen  that 
the  importance  of  fairing  shape  decreases  correspondingly. 

When  blade  thickness  reaches  about  four  percent,  there  is 
almost  no  discernible  difference  due  to  fairing  shape.  To 
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compare  the  relative  effects  of  lairing  shape  and  thickness 
changes,  consider  the  follovring  example »  Select  a  blade  thick 
ness  of  Vl%  with  K  “  .010,  which  is  a  relatively  poor  fairing 
shape  with  a  blxint  nose.  For  this  shape  at  the  proper  solid¬ 
ity,  a  turning  parameter  of  .74  is  possible  without  separation 
If  K  =*  .990  the  shape  is  nearly  the  best  possible. 

Making  this  change  increases  the  permissible  turning 
parameter  to  .S3,  or  a  125^  improvement.  This  same  improvement 
could  be  obtained  with  the  original  shape  by  decreasing  the 
thickness  to  a  little  less  than  and  increasing  the  solidity 
slightly.  If  the  blades  could  be  decreased  to  4^,  with  the 
corresponding  increase  in  solidity,  the  turning  parameter 
could  be  increased  to  .92  without  changing  the  loading.  This 
is  a  24/S  improvement,  or  double  the  best  improvement  that 
could  be  obtained  by  fairing  shape  changes  alone.  As  a  limit¬ 
ing  value,  if  the  blades  could  be  made  infinite sraally  thin, 
the  increase  in  turning  parameter  could  be  as  high  as  50^. 

It  must  be  remembered  that  these  performance  figures  are 
for  cascades  operating  only  at  the  ideal  angle  of  attack.  It 
may  well  be  that  these  increases  in  peak  performance  are  at¬ 
tained  with  some  sacrifices  in  flexibility.  This  would  be  of 
no  consequence  in  a  constant  output  device,  but  for  general 
application,  the  problem  of  operating  off-design  must  be  con¬ 
sidered.  The  analytical  method  described  in  this  study  can 
be  modified  to  consider  the  case  of  operating  at  angles  of  at¬ 
tack  other  than  the  ideal.  This  matter  is  discussed  more 
fully  in  the  section  on  suggestions  for  further  work. 
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5.0  Conclusions 

It  is  concluded  that  the  method  of  generating  cascades 
as  described  in  this  report  is  a  workable  one,  and  provides 
useful  results.  Cascades  are  generated  as  complete  entities, 
which  eliminates  the  need  for  separating  the  effects  of  the 
various  design  parameters.  Use  of  a  high  speed  digital  com¬ 
puter  means  that  large  numbers  of  cascades  can  be  generated 
and  analyzed  quickly.  This  allows  exploration  over  a  wide 
range  of  design  parameters. 

Solidity  in  experimental  cascades  has  often  been  arbi¬ 
trary.  This  investigation  has  revealed  that  the  matter  of 
solidity  is  of  paramount  importance.  For  any  combination  of 
fairing  shape,  thickness  and  camber,  there  is  a  particular 
value  of  solidity  which  will  result  in  the  lowest  loading  of 
the  boundary  layer.  This  loading  can  increase  quite  rapidly 
as  the  solidity  is  varied  from  the  proper  value. 

Data  obtained  in  this  investigation  has  yielded  the  un¬ 
expected  result  that  the  optimum  value  of  solidity  referred 
to  above  is  relatively  independent  of  the  fairing  shape  of  the 
blades.  An  important  application  of  this  result  will  be  in 
planning  experimental  investigations  of  cascade  performance. 

The  actual  values  of  solidity  given  in  this  report  may  be  used, 
with  the  probability  that  they  will  be  close  to  the  optimum 
values  for  any  sort  of  physical  blade  shape. 

Three  parameters  intended  as  indicators  of  the  tendency 
toward  flow  separation  were  compared  in  this  report.  As  a 
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result  of  this  comparison,  it  is  concluded  that  the  parameter 
S,  as  defined  in  this  report,  is  a  useful  indicator  in  this 
respect.  It  is  further  concluded  that  the  value  S  *=  4*25  is 
an  acceptable  design  limit,  which  is  as  consistent  with  NASA 
data  as  this  data  is  consistent  with  itself. 

With  respect  to  the  problem  of  increasing  the  performance 
of  airfoils  in  cascades,  it  is  concluded  that  the  matter  of 
blade  thickness  is  of  much  greater  importance  than  the  fairing 
shape .  The  example  was  given  of  a  cascade  with  blades  of  12% 
thickness.  Possible  increases  in  turning  parameter  for  this 
cascade  due  to  improvements  in  fairing  shape  were  shown  to  be 
on  the  order  of  12%o.  Increases  in  turning  parameter  produced 
by  decreasing  the  blade  thickness  to  4%  with  a  corresponding 
increase  in  solidity,  were  shown  to  be  24%,  or  twice  the  pre¬ 
vious  amount.  The  absolute  limit  of  performance  increase  by 
this  method  would  be  given  by  the  case  of  negligibly  thin 
blaaes.  For  this  case,  the  total  increase  in  turning  parameter 
would  be  50%, 

It  was  pointed  out  in  this  study  that  increasing  the 
performance  of  cascades  by  using  thinner  blades  and  higher 
solidities  would  probably  reduce  the  flexibility  of  performance. 
The  effects  of  operating  off -design  were  not  a  part  of  this 
study  since  the  method  of  generating  the  cascades  produced 
airfoils  operating  only  at  the  ideal  angle  of  attack.  It  was 
stated  that  modifications  could  be  made  in  the  method  to  con¬ 
sider  the  off-design  case,  but  this  was  not  done  due  to  time 
limitations. 
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6.0  Suggestions  for  Future  Work 


This  investigation  represents  only  a  beginning  in  the 
field  of  theoretical  design.  With  the  basic  method  and  com¬ 
puter  programs  available,  it  shoxild  be  a  simple  matter  to  ex¬ 
tend  the  scope  of  investigations  similar  to  this  one.  Inves¬ 
tigations  ought  to  be  conducted  for  a  number  of  average  flow 
angles,  so  as  to  cover  the  whole  range  of  possibilities  in 
cascade  design.  Reference  2  for  example,  has  some  work  on  an 
average  flow  angle  of  zero,  which  represents  the  special  case 
of  an  impulse  turbine.  There  is  no  reason  to  prevent  con¬ 
sidering  negative  average  flow  angles  either.  This  would  be 
the  case  of  turbines  rather  than  compressors,  but  the  same 
method  and  computer  program  would  apply. 

One  might  suppose  that  since  the  flow  is  accelerated 
through  a  turbine,  the  boundary  layer  experiences  a  negative 
pressure  gradient  and  separation  is  not  a  problem.  This  is 
not  exactly  true,  however.  There  are  actually  two  pressure 
differences  to  consider  in  flow  through  a  cascade.  One  is  the 
effect  perpendicular  to  the  stream,  which  is  a  result  of  the 
lift  being  developed  by  each  blade.  The  other  is  the  longi¬ 
tudinal  effect  parallel  to  the  stream  which  is  the  result  of 
differences  in  pressure  from  the  front  to  the  rear  of  the 
blade  row.  Flow  separation  depends  on  the  combination  of  these 
two  effects.  Thus,  in  decelerating  cascades,  the  longitudinal 
effect  reduces  the  allowable  lift  from  each  blade.  In  tur¬ 
bines,  the  longitudinal  effect  increases  the  allowable  lift. 
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If  the  lift  is  increased  enough  to  overcome  the  aid  of  the 
longitudinal  effect,  flow  separation  will  still  occur.  A 
series  of  studies  of  this  matter  could  establish  the  exact 
relation  between  these  two  gradients  and  the  occurrence  of 
separation. 

The  analytical  methods  of  this  report  would  be  a  help  in 
studying  the  occurrence  of  separation  through  solution  of  the 
boundary  layer  equations.  The  detailed  pressure  distributions 
provided  by  the  computer  at  least  open  an  avenue  toward  solu¬ 
tion  of  the  boundary  layer  equations  by  direct  integration. 

In  any  future  work,  it  would  be  convenient  to  consider 
the  velocity  at  more  points  near  the  leading  edge  of  the  blades. 
This  would  allow  more  accurate  determination  of  the  actual 
maximum  velocity,  and  should  eliminate  some  of  the  scatter 
found  in  the  data  for  this  investigation.  The  necessary  modi¬ 
fications  to  the  computer  program  are  not  difficult,  and  only 
lack  of  time  prevented  their  incorporation  in  this  investiga¬ 
tion. 

Another  area  of  investigation  would  be  a  change  in  the 
basic  transformation  which  would  allow  a  finite  angle  at  the 
trailing  edge  of  the  blades.  This  would  make  the  transforma¬ 
tions  similar  to  the  Karmann-Trefftz  transformations  then, 
rather  than  the  Joukowski.  Eliminating  the  cusp  or  radius  at 
the  trailing  edge  then,  would  perhaps  yield  blade  shapes  which 
would  be  more  physically  useful.  As  pointed  out  earlier,  some 
of  the  blade  shapes  developed  in  this  investigation,  while 
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mathematically  useful,  were  physically  poor  because  of  the 
rounded  trailing  edge.  These  improved  transformations  would 
be  particularly  useful  in  the  case  of  symimetrical  impulse 
blading,  where  a  finite  angle  is  desired  at  both  the  leading 
and  the  trailing  edge. 

As  mentioned  previously,  it  would  be  interesting  and  per¬ 
haps  useful  to  consider  the  off-design  case,  where  the  blades 
are  operating  at  some  angle  of  attack  other  than  the  ideal. 

In  principle,  the  method  would  be  to  generate  the  airfoils, 
then  go  back  to  the  original  complex  potential  and  change  the 
average  flow  angle.  Keeping  the  same  transformations  would 
preserve  the  blade  shape,  but  the  circulation  would  have  to 
be  adjusted  to  keep  the  Kutta  condition.  The  idea  seems  simple 
enough,  but  no  doubt  there  will  be  many  difficulties  with  the 
details. 

Still  further  work  could  be  done  in  investigating  other 
shapes  for  the  camber  line  of  the  blades.  Generally,  the 
camber  line  for  the  blades  developed  in  this  investigation  was 
the  slightly  flattened  arc  of  a  circle.  The  generation  of 
other  camber  lines  involves  the  expression  for  the  flow  poten¬ 
tial  in  the  near  circle  plane.  In  this  plane,  a  single  doub¬ 
let  generates  the  flow  contour  which  transforms  into  one  of 
the  blades  of  the  cascade.  If  instead,  a  series  of  doublets 
were  used;  by  varying  the  strength  and  distribution  of  these 
doublets,  great  control  could  be  exercised  over  the  camber 
line  and  blade  shape  of  the  airfoils.  The  ultimate  in  this 
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direction  would  be  where  the  input  is  a  desired  pressure 
distribution,  with  the  proper  camber  and  blade  shape  being 
automatically  produced.  The  resultant  family  of  cascades 
would  be  as  important  a  milestone  to  cascade  design  as  the 
NASA  6-series  airfoils  were  to  isolated  airfoil  design. 
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Boundary  Layer  Loading  Parameter  vs.  Solidity 

For  Constant  Thickness  Airfoils 
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Boundary  Layer  Loading  Parameter  vs.  Solidity 
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Boundary  Layer  Loading  Parameter  vs.  Solidity 
Showing  the  Effect  of  Parameter  K  Variations 
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Design  Solidity 
Jb^igure  13. 


Boundary  Layer  Loading  Parameter  vs.  Design  Solidity 
Bor  Various  Cascades  Beta  =  43°  K  =  .500 


47 


Bo  urtdS  L  C  r  L  o  ^  P<^  r  3  a  t  c  k 


\ 


\ 


r 

'  Gerieo  Tv/o 


9.0 


..•4 


eK" 

>1  8.0 


7.0 

6.0 

5.0 


4.0 


5.0 


2.0 


1.0 


.7  .9  1.1  1.5  1.5  1.7  1.9  2.1  2,5 


Design  Solidity 


j^’igure  14, 

Boundary  Later  Loading  Parameter  vs.  Design  Solidity 
Por  Various  Cascades  Beta  =  45°  a  =  .745 
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Series  Tv/o 
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i’igure  15. 

Boundary  Layer  Loading  laraaieter  vs.  Design  Solidity 
jc'or  Various  Cascades  Beta  =  45°  K  =  .990 
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Series  Two 


figure  16. 

boundary  Layer  Loading  Parameter  vs.  uptimima  Solidity 
for  Various  Cascades  Beta  =  45°  K  =  .010 
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beries  'I'wo 


Optimvun  Solidity 
j?igure  17. 

Boundary  Layer  Loading  Parameter  vs.  Optimvim  Solidity 
Por  Various  Cascades  Beta  =  45°  K  =  .245 
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jb'igure  16. 

boundary  layer  Loading’  iarameter  vs.  Optimum  Solidity 
lor  Various  Cascades  leta  =  45°  K  =  .500 
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Series  Two 
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ligure  19. 

Boundary  Layer  noading  parameter  vs.  Optimum  Solidity 
i'or  Various  Cascades  Beta  =  45*^  K  =  .745 
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3  o  iJ  f)  o'  .J  LaS<'Ji  Pa  t- a  r»  il  t.  c 


Series  Two 


Optimum  Solidity 


i'igure  20. 

Boundary  layer  Loading  Parameter  vs.  Optimum  Solidity 
Por  Various  Cascades  Beta  =  45°  K  =  .990 
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d 'f'  ^  (it  ±  c  ^ 


J^  igure  21. 

Turning  rarameter  vs.  uesign  solidity 

Jj'or  boundary  nayer  Loading  parameter  of  4.25 

(  Corresponds  to  ImaSa  local  Liffusion  factor  of  .5 
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) 


i'igure  22, 

Turning  Parameter  vs  Optimum  Oolidity 

Tor  Boundary  Layer  Loading  parameter  of  4.2b 

(  Corresponds  to  BASn  Local  Liffusion  Tactor  of  .5  ) 
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i-iATHEMTICAL  AI^aLYSIS  uJj'  Thii  i-ic-Triui;  uA  DiiSlGl^ 


I'his  appendix  describes  the  general  mathematical 
developraent  ol  the  method  of  designing  airfoils  in  cascaie. 
Some  of  the  fine  details,  particolarly  on  iterative  proceed- 
ures,  have  been  omitted,  however,  the  entire  calculation 
proceedure  is  given  in  the  program  in  the  next  appen¬ 

dix,  i'his,  along  with  the  explanations  given  in  the  table 
of  symbols,  should  allow  analysis  of  the  method  to  whatever 
extent  as  necessary. 

Basic  xcelations  in  the  hear  Circle  rlane 

ihe  complex  potential  for  the  flow  around  the  infinite 
row  of  closed  bodies  in  the  near  circle  plane  may  be  written 
in  the  form: 

w  =  +  i'i'  =  e  +  (a  +  iB)  Bj^(z)  Ala 

with  the  corresponding  comijlex  velocity; 

dw  =  V  -  iv  =  e~^^  +  (a  +  iB)  Bp(z)  Alb 

dz  ^  ^ 

where 

i'^(z)  =  hj^(x,y)  +  iJ^(x,y)  =  i  coth  |  A2a 

F^(z)  =  R2(x,y)  +  iJ2(x,y)  = 

=  coth  I  J  =  -  i  csch^  I  A2b 

By  expanding  these  functions  and  separating  txiem  into 
their  real  and  imaginary  parts,  it  is  found  tnat; 


b7 


hjL(x,y)  = 

1  /sinh  X  i 

^  ^cosh  X  -  cos  y  j 

A3  Si 

J2_(x,y)  = 

.  ±  I _ §ln,.y _ 1 

^  \  cosh  X  -  cos  y  / 

A3b 

1  (cosn  X  cos  y  -  l) 

~  ~  ^  ^cosh  X  -  cos  y^*" 

kJ)Q 

^  1  (sinh  X  sin  y  ) 

~  ^  fcosh  X  -  cos  y)Z 

A5d 

xhe  expressions  for  velocity  potential,  stream  function, 
and  velocity  components  now  reduce  to: 

=  X  cos/5  +  y  sin +  AKj^(x,y)  -  BJj^(x,y)  A4a 

y  =  -  X  sin/S  +  y  cosfS  +  AJj^(x,y}  +  BRj^(x,y)  A4b 

v^  =  cos/3  +  iut2(x,y)  -  -^^c 

Vy  =  sinjS  -  AJ2(x,y)  -  Bfi2(x,y)  A4d 

i'hese  relations  represent  the  flow  without  circulation 
about  an  infinite  series  of  nearly  circular  bodies  spaced 
at  equal  intervals  S  =  2/r  along  the  y  axis.  I'he  contour  of 
a  typical  body  is  shown  in  iigure  b. 

Location  of  Stagnation  Points 

hue  to  the  interference  caused  by  tho  presence  of  the 
infinite  stack  of  bodies,  the  stagnation  points  on  any  body 
are  displaced  slightly  from  the  position  they  would  have 
occupied  if  the  body  were  isolated  in  the  flow,  r’or  this 
reason,  some  computation  is  necessary  to  locate  these  points. 

^8 


ijesijjnatinji,  t,lu  fz-out  ana  rear  sta£;,riation  ^joints  by  the 
tn.bscripts  a  and  n  respectively,  ana  noting  tlie  type  of  polar 
symriietry  involved,  the  follo^^ring  f;.cts  are  observed,  nt  the 
tv/o  stagnation  points,  both  velocity  components  vanish  by 
aefir.it iun.  aIso,  by  reason  of  the  odd.  symmetry  involved, 

'"a  ^  ~  ~  ^A  ““^B  ‘  However,  since 

the  points  .a  and  B  lie  on  the  same  streamline,  ~  • 

xhe  only  way  in  which  these  two  conditions  on  f'  can  be  satis¬ 
fied  is  if  W,  =  -  0  .  bince  the  stream  function  vanish- 

A  B— 

es  at  the  st^^gnation  points,  it  also  vanishes  at  all  other 

points  on  the  body,  i  ince  tney  are  on  the  s^nne  streamline. 

The  velocity  potential <p  on  the  other  hand  varies  from 

a  maximum  negative  value  at  A  to  a  maxjmum  positive  value  at 

B.  The  value  ^  then  defines  the  amplitude  of  the  func- 

ti  n  on  the  contour  of  the  body.  This  amplitude  may  be 

arbitrarily  prescribed,  and  will  then  fix  the  size  of  the 

near  circle  contour  in  relation  to  the  fixed  spacing  a  =  2T7'. 

Let  im  (x^,yv,)  =  xu  o  so  forth.  Then,  at  the  point 

1  B  B  TxJ 

B  from  me  foregoing  discussion; 


cosfi  +  yB  sinf-i  +  =  arbitrary  value 


^  = 

-Xg  sin^  +  yB  cos/^ 

+  AJ^b  + 

=  0 

A5b 

^x 

cos^ 

+  Alt2B  -  ^'^2B 

=  0 

A5c 

V  = 

y 

sin  i'; 

?  -  Ad2B  ~ 

=  0 

A5d 
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‘i'he  last  tv/o  oi  the  above  equa^j-oiis  may  be  solved  simul¬ 
taneously  for  the  unlcnown  constants  a  and  B  v/ith  the  results; 


A  = 


B  = 


^2B  ^2B 

^2B  ^2B 

+  ^2B  ‘^2B 

^2B  ■*'  ‘^2B 


sin 


cos  p 


A6a 


A6b 


substituting  these  results  into  the  expressions  for  ^  and 
yields  two  equations  in  wnich  the  two  unknowns  are  the 
coordinates  and  y^  of  tne  aft  stagnation  point.  The 
equations  are  trancenaental ,  hovjever,  so  that  an  iteration 
technique  is  necessary  for  solution.  Care  v/as  necessary  in 
this  area  as  all  the  functions  are  multi-valued  in  y,  and  the 
result  could  converge  to  a  point  aiffeiing  by  some  integral 
multiple  of  2  TT  in  the  y  dimension  from  the  intended  point. 

Ihe  general  method  involves  starting  with  low  values  of 
in  effect  close  to  the  isolated  case  vniere  the  location  of 
the  stagnation  points  :s  knov/n,  anu  v/orking  up  to  the  desired 
value  by  modest  increments. 


■Location  of  Points  on  the  Contour 

^nce  the  location  of  tne  aft  stagnation  point  has  been 
found,  the  constants  n  ana  a  can  be  evalu.;.ted  from  equations 
iiSs.  and  Abb  .  xhe  other  points  on  the  contour  can  then  be 
located  from  equations  A5a  and  A5b  by  setting  </>  =  0  and 
assigning  a  suitable  value  to  fp  .  I'nis  computer  program 
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divides  the  interval  from  to  -r  into  20  eq^ual  parts. 

This  yields  20  points  along  each  surface  of  the  airfoil,  or 
40  points  in  all.  As  pointed  out  in  the  text,  it  would  be 
helpful  to  subdivide  some  of  these  intervals  near  the  leading 
edge  to  get  a  better  check  on  the  maximum  velocity  in  this 
critical  region. 

ihe  equations  for  the  location  of  these  contour  points 
are  again  trancendental ,  and  a  similar  iterative  technique 
is  necessary  for  solution,  once  the  coordinates  of  a  point 
are  knovm,  the  velocity  components  are  found  from  equations 
A5c  and  ii.5d.  I'he  final  velocity  v  and  the  slope  c<  of  the 
contour  at  that  point  are  found  from  the  relations; 


A8a 


V 

o(.  =  arctan 

X 


A8b 


At  the  stagnation  points,  this  relation  becomes  indeterminate 
and  L'Hopital's  rule  is  used  for  a  solution. 

basic  relations  in  the  Circle  rlane 

It  is  well  knoxra  that  the  flow  around  an  infinite  rov; 
of  arbitrary  but  identical  contours  can  in  principle  be 
mapped  into  a  corresponding  flow  about  a  single  unit  circle. 
If  the  flov'/  is  of  zero  circulation,  and  the  direction  of  flow 
far  from  the  contours  is  then  the  corresponding  flow  in 
the  circle  plane  is  expressed  bp  tno  following  complex  po¬ 
tential  and  velocity; 
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oi).litting  these  into  real  and  imaginary  parts,  and  placing  r 
=a  =  1  leads  to  the  result  that  along  the  unit  circle, 


V 


K 


=  0 


^  =  constant 


AlOa 

AlOb 


Ihis  verifies  tiiut  the  circle  is  in  fact  a  streamline, 
ij'urther,  along  this  circle: 


cosh(f^  +  cos  ^ 

COShi;?''-  cos^ 


Alla 


V9  -  ^ 


cos 


/^cosh(^  sini^  +  sin  sinh  ^  cos  9 


cosh  2Y'  -  cos  2  & 


Allb 


Singularities  and  stagnation  points  in  the  Circle  Plane 

fo  find  the  stagnation  point  h,  set  =  0  .  I'hen 

equation  Allb  gives: 

tan  =  tan  tanh  )f^  a12 

Substituting  this  into  equation  ^llii  produces  a  result  which 
can  be  reduced  to  the  form: 
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=  cos/1  In  ^  sin/i.2  arctan  Al5a 


wnere 


y  =  y  cosh  2^  +  cos  2/^) 


Alpb 


i’hese  equations  then  must  be  solved  for  T  then  ^  for 
assigned  values  of  ^  and  which  must  be  the  same  as  the 
values  for  and  previously  employed  in  the  near  circle 
plane.  Since  the  functions  are  tranceudental,  an  iterative 
tv-clmique  is  again  required.  Nev;ton's  metnod  is  employed, 


which  involves  calculation  of  tlie  derivative 
noting  this  by  the  symbol  then; 


o)  r 


be- 


-8'i' 


.2 


= _ 

®  ^  2T^  -  cos  2 


A14 


Applying  the  equations  in  this  form  however,  requires  an 
excessive  number  of  significant  figures  in  the  evaluation  of 
T  .  Xo  circumvent  this  problem,  the  following  subltitution 
was  used  to  eliminate  'i'. 


i'  =  cos 


Al5 


I’he  x-’X'evious  equation  then  become 


=  cos  ^ 


2  cos  +/iT 


sin 


sln/?.2  arctanj^ 


^  1 
5  +alj 


Al6a 
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Aiet 


if)'  /2  ]fc03p>  \  (  COB  +^T) _ 

’B  ~(^T/(2  cos/3 +^iy  £1+^2  cos/'««  +  ^'i’)AX7 

Xhe  problem  of  finding  1  for  assigned  values  of  and^ 
b^'  i^ewton's  method  has  no  difficulty  with  regard  to  signifi¬ 
cant  figures,  ‘x'he  parameter^  and  its  hyperbolic  f^.nctions 
may  now  be  found  from: 


sinh  ^  =  j/l'  (  2  cos^  +  /J1)  =  ij (  cosh  2^-1) 


A17 


1‘he  initial  estimate  foi-  newton's  method  is  obtained  by 
neglecting  T  in  the  last  term  of  equation  alba  .  the 
resulting  equation  may  be  solved  explicitly  for  T  . 


v/here 


A  T  = 


2  cos 


1 


2  sin 

cos  /2 


) 


A18a 


A18b 


‘these  expressions  provide  an  excellent  first  estimate  for 
1‘,  and  in  fact  for  the  special  case  Qt^=-  0  ,  they  are 
exact. 

with  ^  and  ^  known,  the  location  of  any  point  <9. 
having  an  assigned  value  of  is  defined  implicitly  by 

equation  ^illa.  I'o  establish  tiie  location  of  points  corres¬ 
ponding  to  those  in  the  near  circle  plane,  all  that  is  neces¬ 
sary  is  to  have  ti.e  respective  values  of  be;  the  same  in 
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the  ti70  planes,  as  before,  an  iterative  solution  is  re¬ 
quired  . 

Circulation 

'ihe  complex  potential  ana  velocity  for  a  purely  cir¬ 
culatory  flow  of  unit  circulation  in  the  circle  plane  are 
as  follows: 


Aw  =  - 


i  cos 

2 


In 


A19a 


^w)  ^  -  2i  cos  sinh  ^  S 

-  2  oosh/5.  j2+l 


A19b 


i  ^5? 


r  e 


A19c 


Separating  these  into  real  and  imaginary  parts,  it  is  again 
found  that  the  unit  circle  r  =  1  is  a  streamline,  along 
whiph  the  velocity  distribution  is  given  by 

.  Y  _  -  s  inh  2  ^ _  a20 

^  cosh  2^  -  cos  2 


This  circulatory  component  must  be  evaluated  at  the  same 
points  previously  establishec^round  the  periphery  of  the 
circle.  By  superposition,  the  total  velocity  at  these  points 
is  then  given  by; 


+  jj 

0  ^ 
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otaf:nativn  j.-‘oints  with  circulatiuri  ii>  Circle  xlaiie 


Just  as  the  stagnation  points  ana  B  are  defined  by 
the  condition  V  =  0  ,  so  are  the  stagnation  points  L  and  9? 
defined  the  corresponding  condition  0  .  introducing 

this  condition  into  the  foregoing  equations  yields  the 
essential  result 

-  (5>„  =  S’,  -  (9,  =  r  =  aroslnr?a  . 

®  ^  ^  cosh^  oos  (i 

A22 

urice  the  points  B  and  'i  have  been  located  in  the  circle 
plane,  the  corresponding  values  of  and  may  be  found 
from  equation  Alla,  introducing  these  \alues  into  the  near 
cii’cle  plane  leads  to  the  location  of  tiie  points  L  and  'i  in 
that  plane  by  solving  equations  A4a  and  A4b. 

manoinif^'  -uetyeen  Circle  and  hear  Circle  Planes 

Happing  between  the  circle  cr  J  plane  and  the  near 
circle  or  a  plane,  as  implied  in  tne  preceeding  discussion, 
is  accomplished  despite  the  fact  that  no  explicit  function 
for  a  in  terms  of  ^  or  vice-versa  is  Icnovm.  -^he  mapping  is 
dene  implicitly  through  the  fact  that  the  non-circulatory 
flows  in  the  two  planes  must  correspond,  or  in  other  v^ords, 
tiirough  the  condition  that 

wi,  s)  =  w  t  5 ) 


A2pa 


i^ifl'erentiatini^  this 


dw  I 
djs| 


But 


I dw  I  =  V 

I  dz  1 


A25b 


-‘23  c 


and 


A23d 


Hence  the  mapping  derivative  becomes 


xhe  flew  with  circulation  has  been  completely  defined 
in  the  foregoing  development  for  tile  circle  plane.  In  the 
near  circle  plane,  however,  only  the  non-circulatory  flow 
has  been  explicitly  defined.  explained  in  the  text,  add¬ 

ing  circulation  terms  to  the  complex  potential  in  the  near 
circle  plane  produces  unv/anted  distortions.  Therefore,  the 
method  used  is  to  map  the  known  flow  v^rith  circulation  in  the 
circle  plane  back  into  the  near  circle  plane  using  the  map¬ 
ping  relations  previously  established.  These  mapping  rela¬ 
tions  are  of  course  independent  of  circulation,  •‘•hus,  for 
the  final  flow  v/ith  circulation 

w^  (z)  =  (5)  =  ^3)  +  B^Aw  (5)  A24a 
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itifferGutiating  this  gives 


II 

o 

\n* 

B 

Q  d  / 

A2^b 

dz  I 

1—  1 

|d5( 


Substituting  in  terms  of  velocities,  and.  introducing  proper 
algebraic  signs  gives  finally  for  the  velocity  in  the  near 
circle  plane; 


V  =  yv. 

V 

V 


=  V  +  B 

o 


A25 


Basic  helations  in  the  Cascade  Plane 

The  basic  idea  of  this  entire  method  is  to  map  the  near 

1 

circles  into  a  cascade  of  airfoils  in  much  the  same  manner 
that  the  doukov;’ski  transformation  may  be  used  to  map  a  single 
circle  into  an  airfoil.  The  equations  for  the  transformation 
are  identical  in  form  to  those  already  developed  for  the 
flow  in  the  near  circle  plane,  however,  all  coordinates  are 
evaluated  with  respect  to  axes  shifted  slightly  from  the 
original  position,  denoting  these  shifted  coordinates  by 
primes,  the  maxjping  function  employed  is  specifically; 

* 

=  e"^^z'  +  (A'  +  iB')  A25a 

-4^  =  u  -  iu  =  e"^^'  +  (A'  +  iB')  Fp(z')  A25b 

dz '  X  y 
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where 


' )  =  T  coth  |-- 
F^iz')  =  -i  csch^  ~ 

hote  that  the  above  functions  are  exactly  of  tne  same  form 
as  equations  A1  and  A2  .  All  subsequent  results  may  readily 
be  obtained  through  this  analogy. 

Determination  of  Airfoil  ohape 

In  order  to  transform  tnc  near  circle  into  the  required 
airfoil  shape,  it  is  necessary  to  locate  the  singular  points 
i  and  ^  of  the  mapping  transformation  in  the  proper  relation 
to  the  stagnation  points  a  and  r  of  the  flow,  normally,  r 
is  locatea  a  short  distance  inside  the  contour,  just  opposite 
L  on  the  surface.  Similarly,  point  lies  inside  the  con¬ 
tour  on  a  line  normal  to  the  contour  at  T.  nefer  to  i'igure 

5,  rhe  significance  of  ti^ese  locations  may  be  qualitatively 

PL 

understood  from  the  fact  that  as  the  relative  distance  — 

=  ^  is  made  small,  the  radius  of  curvature  at  the  nose  of 
the  airfoil  decreases  rapidly,  and  iii  the  limit,  as  ^  =  0 
tne  nose  becomes  cusped.  Similarly,  the  relative  distance 
=  pK  controls  tne  tail  curvature.  In  the  limit  as 
=  0  the  tail  becomes  cusped  in  the  manner  of  a  uoukow- 
ski  airfoil,  fhe  origin  of  coordinates  for  the  transforma¬ 
tion  io  shifted  from  the  origiiial  location,  the  midpoint 


A26a 

A26b 
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between  points  a  and.  13  ,  to  the  nev;  location,  the  midpoint 
between  points  P  and  Q, 


Velocity  and  i^lope  ^t'long  the  Airfoil 

The  complex  potential  including  circulation  ii^  the  cas¬ 
cade  and  near  circle  planes  may  be  equated 

w  {t*)  =  w  Iz)  A27a 


differentiating 


net 


conjugate  of  velocity 
on  airfoil 


conjugate  of  velocity 
on  near  circle 


=  1 


A27b 


A27c 


A27d 


A27e 

..27f 


Prom  which,  ax  ter  separating  amplitudes  and  phase  angles; 


V  = 


V 


(h) 

u 


A28a 


A28b 
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Fortran  Program  "Cascade"  and  Operating  instructions 

Jn  the  following  pages  is  tne  i'unliuu'f  program  "Cascade" 
used  in  deriving  the  data  used  in  this  study,  input  data  is 
by  means  of  data  cards  inserted  at  the  end  of  the  program, 

*iny  nwriber  of  data  cards  may  be  used.  J-he  first  card  indica¬ 
tes  how  many  additional  data  cards  are  to  be  read. 

bach  additional  card  lists  the  number  of  values  of  the 
parameters  ,  ,  ,  and  K  to  be  computed,  their 

initial  values,  and  the  increment  by  o'hich  each  value  increa¬ 
ses.  also  listed  are  the  print  control  quantities  hux^Jti  and 
i'ihAihiL,  which  control  the  extent  of  the  print  out  as  indica¬ 
ted  in  the  table  of  symbols,  '-i-he  computei  then  calculates 
results  for  every  possible  combination  of  the  five  parameters, 
prints  the  results,  and  then  proceeds  to  the  next  card. 

'ihe  standard  print  out  contains  the  five  inpvt  parameters, 
the  solidity  and  stagger  angle  of  the  cascade,  ihe  airfoil 
thickness,  and  the  coordinates,  slope,  velocity,  and  boundary 
layer  parameter  for  each  of  the  42  points  of  Lhe  airfoil. 

Condensed  Print  uut 

when  surveying  large  numbers  of  cascades,  as  was  acne 
for  this  study,  it  is  convenient  to  condense  the  data  to 
reduce  the  volume  of  output.  Pago  85  shows  the  modifications 
made  to  the  basic  ^^rogram  to  accomplish  this,  briefly,  all 
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of  the  standard  print  conunands  weio  removed.  j.hen  conimands 
to  print  the  heading  only  were  inserted  early  in  the  program 
right  after  changes  in  sigma,  'ihe  modified  print  out  then 
contains  no  airfoil  data  at  all  except  for  the  maximum  value 
of  the  boundary  layer  iD.irameter.  Operation  of  the  program 
is  identical  v;ith  the  basic  program. 

Iterative  Qcheme  for  Constant  -t-hickness  Airfoils 

Page  84  shows  further  modifications  made  to  the  basic 
program;  to  produce  cascades  v;ith  the  desired  airfoil  tnick- 
ness.  In  effect  this  eliminates  the  parameter  ^  since  it 
is  varied  as  necessary  to  produce  the  required  thicknesses. 
Ihe  initial  value  of  £  is  still  read  from  the  data  card,  but 
the  column  which  indicates  the  number  of  £  to  be  calculated 
no  longer  has  significance,  ihe  program  as  written  will 
demand  thicknesses  from  ^^/o  to  24/^0  in  increments,  other 
values  for  thickness  can  be  obtained  by  modifying  the  appro¬ 
priate  commands.  Ko  means  of  doing  this  by  data  cards  is 
provided,  fhe  print  out  is  of  the  condensed  form  described 
previously . 
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,ETl(42}j.  ePI(42),  VVI  (42)  tVTl(42) 
,0DI(12), (TITLE  ( 10) 


OELI(42)fTEI 
,YNI(4 


CIRCLE 


•  BEJ 


1?141593 


•JOB  JENISTA 

PROGRAM  CASCAOe 

OIMENSIQN^XI I42),YI(42),RII42),VIJL2) 

1(42),  DVI(42),  XII(42), - - 

22),XNI(  42),  V2I(42) 

I  D*0 
201  READ  204. IB 

204  FORMAT  (112) 

?,S§iNBE,NSG,NBZ,NEP,NAK,BE1  ,SG1,B21,EP1  ,AK1,BE1,  SGI,  BZI  ,EBI, 

iAr>lf  NUPKf  NcXPR 

205  FORMAT  (512,  10F5. 2,212) 

PART  I A  constants  FOR  NEAR 

BEJ«eEl 
002  J«ltNBE 
BE»  (3.141593/180.) 

SNBE«SINF(BE) 

CNBE*COSF(BE ) 

SGK>SG1 
003  K«1,N 
PHEB«SGK* 

C22>:SNBE 
C23*CNBE 

IF(K>1)  401,401,402 
DRUP-1.5708 
DALP«0. 

SIP»0. 

RUP*0. 

ALP*BE 

IF( SGK-SIP-.5)  404,404,403 
SIG*SIP^.5 
RU=RUP+DRUP»  .5 
ALU»ALP+DALP».5 
IF(NEXPR)  S2C,520,518 

WRITE  OUTPUT  TAPE  2, 5 19, ORUP, DALP, SI P, SIG.RUP ,RU , ALP . ALU 
FORMAT!/  7H  DRUP  *F8.4,  5X  6HDALP  «  F9.4,  /  6H  SIP  -F  9.4,  5X 

1  SHSIG  «F10.4,  SX  SHRUP  >F10.4,  5X  4HRU  «F11.4,  SX  SHALP  -MO. 4, 

2  5X  5HALU  -Flo. 4/  ) 

520  CALL  BAKER  ( S IG, XU, YU,RU . ALU, AU,BU, R lU, AJ 1U, R2U, AJ2U , « 00000 1 , 

1  A2U, ITU,C22,C23,NEXPR) 

DRUP»(RU-RUP )/.5 
DALP-(ALU-ALP)/.5 
SIP-SIG 
RUP-RU 
ALP-ALU 

IF(NEXPR)  522,522,521 

WRITE  OUTPUT  TAPE  2, 5 19, ORUP.DALP. SI P, SI G, RUP, RU, ALP , ALU 
CONTINUE 
GO  TO  402 

'RB-RUP^DRUP»(SGK-SIP) 

AL8»ALP4DALP*(SGK-SIP) 

IF(NEXPR)  525,525,523 

WRITE  OUTPUT  TAPE  2, 524, DRUP ,DALP,SI P, SGK, RUP, RB, ALP , ALB 
FORMAT!/  7H  DRUP  -F8.4,  5X  6HDALP  -  F9.4,  /  6H  SIP  *F  9.4,  5X 
1  5HSGK^«F10 .4^  5X  5HRUP  -F10.4,  5X  4HRB  -F11.4,  5X  5HALP  -F10.4, 

525^CALL  baker  (*"  °*SGK, XB,YB,RB, ALB, A, B, R1B,AJ1  B,R2B,AJ2B,. 000001, 

1  A2  -  - 


401 


402 

403 


518 

519 


521 

522 

4  04 


523 

524 


.  BAKER  (  SGK,XB,) 
k2,IT,C22,C23,NEXPR) 
>«(RB-RUP )/( SGK-SIP) 


526 

527 


DRUP»  .  _  _ 

DALP- (ALB-ALP) /(SGK-SIP) 

SIP-SGK 
RUP-RB 

ALP-ALB  75756 

WRI Ti^OUTPUT^f APE* 15524, DRUP, DALP, SI P, SGK, RUP, RB, ALP ,ALB 
CONTINUE 
SNALB-YB/RB 
CNALB-XB/RB  ,  .  „ 

ALBO  «( 180. /3. 141593)* 

XA— XB 
YA— YB 


ALB 
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ALA«ALB+3. 141593 

RA*RB 

XI(20)«XA 

XI(40)«XB 

Yn20)»YA 

YI  ( 40)*YB 

RI  (20)>RA 

RI(40)»RB 

ALI ( 20)*ALA 

ALI  (40)=ALB<»6. 283186 

Vn2G)sO. 

VI(40)>0. 

R3B«-2.* (RlB«R2e-AJlB«AJ2B) 

AJ3B«-2.« (AJ16«R2B«rTb*A J2B) 

DVXB*A«R3B-fl*A43B 

0VYB»B*R3B+A«AJ3B 

0V2B*SQRTF(0VXB«0VXB4-0VYB«0VYB) 

AHI (40)*ASINF(DVYB/DVZB) -1.5707965 

AHI (20)«AHI(40)«3.141593 

A22=A«A+B«B 

A23»R1*RUAJ1*AJ1 

A24*R2«R2  4^AJ2«AJ2 

W2ZB»2.*SQRTF  (A22*A23*A24) 

IF  (NOPR) 120.120, 121 

120  WRITE  OUTPUT  TAPE  2,  122,  BE J.SKBE.CNBE .SGK, PHE B, A.B .XB ,Y B, RB, ALBC 

1 .  A2 .  IT 

1220F0RMAT  (IHI///  46X  28HRESULTS  IN  NEAR  CIRCLE  PLANE  // 

1  lOX  6HBETA  »F6.2,  5H  DEG.  8X  lOHSIN  BETA  «  F7.4,  8X 

2  lOHCOS  BETA  »  F7.4,  8X  7HSIGMA  »  F7.3,  / 

3  lOX  BHPHKB)  »  F8.4,  9X  3HA  *  F14.4,  8X  3HB  «  F14.4,  // 

4  12X  IHI  9X  IHX  IIX  IHY  llX  1HR  5X  12HLAHBDA,  DEG.  6X  1 HV 

5  6X  IIHALPHA.  DEG.  3X5HERR0R  5X  lOHiTS  SHIFT  // 

6  12X  IHO  3F12.4,  F12.2,  6X  6HO.OOOO  8X4H*  E12.2,  16, 

7  6H  •) 

121  CONTINUE 

PART  IB  POINTS  ON  NEAR  CIRCLE 
DC 


13 


14 


IS 


49 

50 


36 

37 


DC  12  I  » 1,19 
IF(l-l)  13,13,14 
PHN=.9«PHE8 
AL>AL6«.45 
R»RB 

GO  TO  15 
AI»I 

PHN=PHEB«{ I.-AI/IO.) 

AL»AL+( 3. 141593+ALB-AL)/ (21 .-AI ) 

X=R«COSF(AL) 

CALL^SSRixiY.VX, VY,V,C22 ,C23, A, B, PHN, PHEB , ITS, SNALB, CNALB , 1 . , 
1  ITSIDE.ERS.KEXPR) 

XI( I )=X 
YI { I )=Y 

R»SGRTF(X*X4-Y«Y) 

RI(I)»R 
V I  ( I ) =V 

AHI { I )»ATANF(VY/VX) 

IF(VX)49,50,50 

AHK  I)»AHI(I)+3.141593 

CONTINUE 

AHO  »( 180. /3. 141593)*  AHI ( I ) 

AL*ATANF(Y/X ) 

IF( X)36,37,37 
AL»AL*3. 141593 
CONTINUE 

ALO  «( 180. /3. 141593)*  AL 

ALI ( I )»AL 

XI<I*20)»-(XI(I)) 

YI(  I+20)»-(YI (I ) ) 

VI(I*20)*-(Vni)),  ,, 

ALK  1*20)* ALI(  I  )*3.  14159  3 
AHI ( 1*20)  «  AHI ( I ) 
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124 

125 
12 


RI  (  U20)  «  Rid) 

IF  tN0PR)_124,124,12 


11-  imurrv#  l<««f  1^ 

i!!SJlJ.?,Oyi?yT,lAPE,  2,  125.  I.X.Y.R.ALn.V.AHO.ERStl 
(n3.3F12.4,  F12.2.M2.4,Fl2v2,El2.2,2li) 


ITS.ITSIOE 


528 

529 

530 


16 


CONTINUE 
PART  IIA  CONSTANTS  FOR  CIRCLE 
IF(NEXPR)  530,530,528 
WRITE  OUTPUT  TAPE  2,  529 
FORMAT  (///) 

IT3»0 

POWER  ■  (PHEe-2.*BE*SN8E )/CNBE 
DET«(2,*CNBE)/(EXPF( POWER)-!.) 

ARC  1  •  1.4:(2.*CNBE)/DET 
ARG2  •  SNBE/ (CNBE^DET) 

PHE«  CN6E*LOGF( ARG1)+SNBE«2.*ATANF(ARG2) 

ER3=PHE-PHEB 
RER3*ABSF( ER3)/PHEB 
IT3*IT3>1 
T=CNBE+0ET 
Q1*CNBE+T 
Q2*l .+Q1«0ET 

PHP*(-2./DET)*{T/Ql)«(T/Q2) 

DETNU  «  DET-ER3  /PHP 
IF(NEXPR)  516,516,505 

WRITE  OUTPUT  TAPE  2,506, IT3, GET, PHE,RER3, DETNU 

FORMAT(6H  IT3  »I4.5X5H0ET  «E10.4,5X5HPHE  «F1 0.6, 5X6HRER3  «E10.4, 
5X7HDETNU  »E1C.4) 


505 

506 


1 


516  IF(.000001-RER3)  17,19,19 

17  IF(IT3-25)  18,19,19 

18  IF(DETNU)  60,61,61 

60  DET=DET/2. 

GO  TO  16 

61  DET*OETNU 
GO  TO  16 

19  CH2G=1.+2.«01*0ET 

SH2G=SQRTF  ( CH2G*CH2G-1 • ) 
SHG  »  SQRTF  ( (CH2G-1 . )/2. ) 
CHG  =  SQRTF  ( (CH2G+1 . ) /2 • ) 
C2=CH2G 
C3=SH2G 
C4*SHG 
C  5  *0  HG 

TEB»ATANFdC4/C5)*(C22/C23) ) 
TEBD  »(180./3. 141593)* 
CNTE8=C0SF(TE6) 

TEA=TEB+3. 141593 
TEI(20)*TEA 
TEI (40)=TEB+6. 283186 
>VTB=0. 

VTA=0. 

0VB=-(C3)/(C2-C0SF{2.*TEB) ) 

0VA*DVB 

VTI(20)=0. 

VTI(40)=0. 

DVI (20)=DVA 

DVI (40)*DVB 

A22»C23*C5*CCSF{TEB) 

A23=C22*C4*SINF(TEB) 

A24«C2-COSF(2.«TEB) 

W2TB=4.*( A22+A23)/A24 
IF  INOPR)  126,126,128 


TEB 


1260WRITE  OUTPUT  fAPE*2,  1 27 , SHG.CHG ,DET , PHPfTEBO.DVB ,RER3» IT 3 
1270FaRHATt//41X37HCORRESPONOING  RESULTS^IN  tlRCLE  PLANE// 


V  THETA 
F20.4. 


10X 


lOX  12HSINH  GAMMA  «  F8.4,  5X  12HCdSH  GAMMA  *  F8.4,  5X  7H0EL  T  ■ 

I?  t'hETA  lOX  13H0ELTA 

5HERR0R  nx  3H1TS  //  1  5X  IHO  F17.2,  15X  6HO.OOOO 

E19.2,  113) 

128  CONTINUE  ^ 

PART  IIB  POINTS  ON  CIRCLE 


1 

2 

3 

4 

5 


75 


r 


DO  20  1=1,19 
AI*I 

PHN=PHEB» ( I.-AI/IO. » 

IF( 1-1)21 ,21 ,22 

21  TEOLO  »  TEB 

♦  SCRTF  (.2»PHEB/W2TB) 

GO  TO  23 

22  TEOLD  »  TE 

TENU  1  »  TE-.UPMEB/VTE 
,  IF  (TEA-TENU  1)  53,53,54 

53  TE  =  (TE<^TEA)/2. 

GO  TO  23 

54  TE  *  TENU  1 

23  IT4-0 

25  C6»SINF(TE) 

C7»C0SF(TE) 

C8*C0SF(2.»TE) 

A2*L0GF((C5*C7)/(C5-C7)) 

A3«2.»ATANF(C6/C4) 

PH=C23»A2+C22»A3 

A2=-(C23*C5»C6)*C22»C4*:7 

VTE=4.*A2/(C2-C8) 

DPH  =  PH-PHN 
IT4*IT4*1 
A2=ABSF (DPH) /PHEB 
IF(N£XPR)  511,511,512 

WRITE  OUTPUT  TAPE  2,  513,  I,  IT4,  TE , DPH, VTE . A2 
FORMAT  (  6H  I  =14,  5X5HIT4  =14,  6X4HTE  «El0.4, 


24 

26 

51 

55 

56 

57 


SX4HA2 
111 
26,  1 


111 


TE 

GO 

IF 

TE 

GO 

TE 

GO 

TED 


51,51  ,55 


512 

513 

1  5X5HVTE  =£10.4, 

511  IF(IT4-25)  24,111, 

IF  (  .000001  -  A2) 

TENU  *  TE-DPK/VTE 
IF  (TENU  -  TEOLD  ) 

-  (TE+TE0L0)/2. 

TO  25 

(TEA-TENU)  56,56,57 
=  (TE*TEA)/2. 

TO  25 
=  TENU 
TO  25 

(180. /3. 141593)* 


cElO.4) 
11,  111 


6XSH0PH  «E10.4, 


TE 


TEI (  I)  =  TE 
VTI (  I  )=VTE 
DVK  I)=-C3/(C2-C8) 

TEI  (  1*20)  *TEI(  I  )*3.1  41593 
VTI ( I*20)»-(VTI  (I  )  ) 

DVI ( 1*20)  «  DVI ( I ) 

IF  (NOPR)  129,129,20 
WRITE  OUTPUT  TAPE  2,  130, I, T( 
FORMAT  (116,  F17.2,  F21.4,  FJ 
CONTINUE 

PARTIIIA  CONSTANTS  FOR  AIRFOIL 
6ZL=BZ1 
D029L=1,NBZ 
BAM=A6SF(BZL) 

IF( BAK-.OOl ) E5,85,86 
TAU=0. 

XT  =  XB 
YT=YB 
XI{42)=XT 
YI(42)=YT 
VT»0. 

AHT=AHI(40) 

AHL=AHI (20) 

AHI (41 )=AHL 
AHI (42 )*AHT 
TET=TEB*6. 283186 
PHT*PHE8 
XL*XA 
YL«YA 


129 

130 
20 


85 


0,VTE,0VI( I) 
:0.4,  E19.2, 


,A2,IT4 
1  13) 


76 


XI(U1)=XL 

YI(UU*YL 

VL»0. 

TEL=TEA 
PHL»-PHEB 
GO  TO  87 
86  CONTINUE 
C2=CH2G 
C3-SH2G 
C4*SHG 
C5=ChG 
C22*SN8e 
C23-CNB6 

TAy«ASINF( (B2L/U.)»( SH2G/CHG)*(CNTE8/CN0En 

TeT»TEB-TAU+6.283186 

TEL*TEA+TAU 

A2=SINF(TET) 

A3*COSF(TET) 

A4=(C5+A3)/(C5-A3) 

A5aA2/CU 

PHT=C23»LOGF(A4) +022*2. •ATANFIAS) 

A2*SINF(TEL) 

A3»C0SF(TEL) 

A4a(C5+A3)/(C5-A3) 

A5=A2/C4 

PHL=C23*LOGF(A4 )+C 22*2 .* ATANF ( A5 ) 

ENTx30.+10.*PHT/PHEB 
1*39 
AI»I 

IF(ENT-AI)  3C.31,31 
1*1-1 
GO  TO  32 

A2»( TET-TEI ( I ) 1/ (TEI (I+l  )-TEI(  I  1 ) 

RT*RI(I)  +  A2*{RI(I  +  1)-RI(  I)) 

ALTxALl  (n  +  A2*(ALI(  I+l  )-ALI(  I  11 
XT»RT*COSF( ALT) 

YT*RT*SINF(ALT) 

ENL*30.+10.*PHL/PHEB 
1*21 
AI«I 

IF(AI-ENL)34,35,35 
1*1+1 
G0T033 

A2*(TEL-TEI( I ) ) / ( TE I ( I+l )-TEI ( I ) ) 

RL*RI(I )  +  A2*(RI(I  +  l)-RI(  I)) 

ALL»ALI (I)+A2*(ALI ( I+l )- ALI ( I ) ) 

XL=RL*COSF( ALL) 

'YL=RL*SINF(ALL ) 

CALL  SSR(XT,YT,VX,VY,VTtC22tC23tA,B,PHTtPHEB, ITStSNALBtCNALBt-1 . t 
1  ITSIDE.ERST.NEXPR) 

XI (42)*XT 
YK  42)*  YT 
AHT*ATANF(VY/VX) 

IF( VX)63,64,«4 
63  AHT*AIIT+3. 141593 

***  CALLSSRUlIyLi  VX,  VY,VLiC22tC23t  AtBtPHLtPHEBt  ITS  tSNALBtCNALBt- 1  .t 


32 

30 

31 


33 

34 

35 


1  ITS 
XI  (41 
YI  (4  1 


65 

66 
87 


IF(VXio;}lootco 

AHL*AHL+3. 14 1593 

AHI(41)xAHL 

CONTINUE 

TNBE*SNBE/CNBE 

TN8E2»TNBE-BZL/2. 

CV2-(l.+TNBE*TNBE)/( 1 .+7 NBE2*TNBE2 ) 
CV«SQRTF(CV2) 

EPH*EP1 


DE,£RSL»NEXPR) 
XL 
YL 


AHL=ATANF(VY/VX) 


65i66»66 


77 


I 


DO  40  M=1,NEP 
AKN=AK1 
DC  4  1  Nal ,NAK 
IF(EPM)  99,98,99 

98  IF(N-n  99,99,41 

99  CONTINUE 
XP=XL-EPM*R0»SINF( AHL) 
YP»YL+EPM*R8«CQSF(AHL) 
XC*XT-AKN*EPP*R8*SINF( AHT ) 
YQ=YT+AKN*EPM*RB*C0SFCAHT) 

XZ*( XP+XQ) /2. 

YZ*(YP+Y0)/2. 

XQP=XQ-XZ 

YQpxYC-YZ 
X=XQP 
Y=YGP 
C2=S INH(X) 

C3=C0SH(X) 

C6=S INF (Y) 

C7=C0SF(Y) 

R1=(C2/(C3-C7) )/2. 
AJ1=(C6/(C3-C7) )/(-2.) 

PX=C3*C7-1 , 

PY=C2*C6 
PZ*PX*PX+PY*PY 
R2=PX/( PZ*(-2. ) ) 

AJ2=PY/ (PZ*2 . ) 

QZ=R2*R2+AJ2«AJ2 

PU*Y*QZ-R2»AJUAJ2*R1 

PV=X»QZ-AJ2»AJ1-R2*R1 

BEP=ATANF(PU/PV) 

C22=SINF(BEP) 

C23=COSF(BEP) 

QX=AJ2*C22-R2*C23 

0Y=R2*C22+AJ2»C23 

AP*OX/QZ 

BP=QY/QZ 

AP22=AP*AP+BP»BP 

AP23=R1»R1+AJ1*AJ1 

AP24=R2*R2+AJ2*AJ2 

WP2ZB=2.*SCRTF(AP22*AP23*AP24) 

PARTIII8  POINTS  ON  AIRFOIL 
1  =  1 

42  X=XI ( I )-XZ 
Y=YI ( I )-YZ 
C2=SINH( X) 

C3=CCSH(X) 

C6=SINF(Y) 

X7=C0SF(  Y) 

Rl=C2/( 2.* (C3-C7)  ) 
AJl=C6/(-2.*(C3-C7)) 

PX=C3*C7-1. 

PY=C2*C6 
PZ=PX»PX+PY*PY 
R2*PX/( PZ* (-2.  )  ) 

AJ2=PY/(2.*PZ)  ,  „  .  , 

XI I  (  I  )*X*C23-»Y»C22  +  AP*R1-BP«AJ1 
ETI ( I )=-X*C22+Y*C23+AP*AJl48P*Rl 
A2=C23+AP*R2-BP*AJ2 
A3=C22-AP*AJ2-0P*R2 
A4=SQRTF ( A2*A2+A3*A3) 
IF(A4-.0001)  171,171,172 
171  RP3=-2.*(R1»R2-AJ1*AJ2) 
AJP3=-2,*( AJ 1*R2+R1*AJ2) 
0UXP=AP*RP3-eP*AJP3 
DUYP=BP*RP3+AP*AJP3 
0UP=SQRTF(0UXP*CUXP+0UYP*0UYP)  ^ 
EPI ( I)=ASINF(OUYP/OyP)-l .5707965 
IF(ALI ( I)-ALI (30)  )  173,173,113 
173  ePI(I)=EPI(I)+3, 141593 


172 

600 
1 13 


88 

89 

38 

39 

43 

44 


59 

62 

46 

47 


GO  TQ  113 

EPI {  I  )  =  ASINF  { A3/A4) 

IF( A2)  600,  1  13,1  13 
EPI{I)=3.1U1593-EPI{I) 

DELI  {n  =  (180./3.141593)«(AHI(n-EPI(  1)1 
A6=ABSF(VTI{  I ) ) 

IF{A6-.0001)  88,88,89 
A5»SQRTF{W2ZE/W2T8) 

GO  TQ  38 
A5«V  n I )/VTI  ( I  ) 

IF(ABSF{A4)-.000l )  39,39,59 
IFIBAM-.OOOI )  43,43,44 
yyl ( I)«SQRTF{W22B/WP22B) 

GO  to  62 

R3*-2.*{R1*R2-'AJ1*AJ2) 

AJ3=-2.*{AJ1  •R2  +  RUAJ2) 

WP2Z=SQRTF{ AP22»{R3*R3+A J3*AJ3) 1 
SHX=SINH( XI{ I ) ) 

CHX=CQSH{XI { I) ) 

SNY=SINF(YI { I) » 

CNY=COSF{Yn  I)  ) 

Rl-( SHX/2. )/(CHX-CNY) 

AJl={-SNY/2. )/{CHX-CNV) 

PX=CHX*CNY-1 . 

PY=SHX*SNY 

PZ*PX*PX+PY*PY 

R2=PX/{-2.*PZ) 

AJ2=PY/{2.*PZ) 

R3— 2.  *  (Rl  *R2-AJ  1*AJ2) 
AJ3=-2.*{AJ1*R2+R1*AJ2) 

W2Z=SQRTF{ {A*A+B*B»*{R3»R3+AJ3*AJ3) ) 
SNTE  =  SINF{TEI(n  ) 

CNTE»COSF (TE  n I ) ) 

SN2TE=SINF{2.*TEI (I) ) 

CN2TE=C0SF{2.*TEI( I ) ) 

B22=CNBE*CHG*CNTE 

B23=SNBE*SHG«SNTE 

B24=CH2G-CN2TE 

W2T=-4.*{ B22+B23)/B24 

DW2T*{2.*SH2G*SN2TE)/(B24«B24» 

A25  =  W2Z-'A5*A5*W2T 
A26={DVn I )»A25)/A6-A5*A5«CW2T 
VVI { I  )  =  {W2Z+BZL*A26)/WP2Z 
GO  TO  62 

VVI{I)  =  {Vni  )+BZL*OVI  (  I)  *A5)/A4 
IF{ I-U2  )  46, 47,47 
1  =  1  +  1 
GQT042 
CONTINUE  ' 

D£LX=XII(42)-XII (41) 

DELY=ETI {U2)-ETI (41 ) 

CHQRD=SCRTF{DELX*OELX+OELY*DELY) 

SQLID=CH0RD/6.283186 

SNA=DELY/CHQRD 

CNA=CELX/CHOPD 

ARQT*( 180./3.141593)*ASINF  (SNA) 
STAG={180./3.141593)*{BEP+ASINF(SNA) ) 


710 

716 

715 


AHTO  ={ 18C./3. 141593) *  AHT 
TETD  =(  18C./3. 141593) *  TET 

AHLD  ={ 1 8C./3. 141593) *  AHL 

TELD  ={ 1 8C./3. 141593) *  TEL 

TAUD  ={ 180. /3. 141593) *  TAU 

BEPD  ={  180. /3. 141593) •  BEP 

IAA=  20 

00  717  ID0»1,  1C 
lAA  =IAA-'l 

t5o=TANF  (DELI{IAA)*3. 141593/180.) 

AA-  XNnn-XNKIAA)  -  (Y  NK  IAA)-YNI  ( 1)  )  *700 
IF  (AA)  714,714,713 
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F6.3t lOXeHB  ZERO 
«  F6.3,// 


713  1  =  1-1 

GO  TC  715 

71U  B8E=ATANF  ((YNKI  +  l)  -YN  I  (  I  )  ) /(  XN  K  I  ^  1  )-XNI  (  I  )  M 
DDEL  =  .5*(BeE+0ELI(IAA)*3.141593/ia0.) 

,,,  CC  =  ABSFl  (YNI(IAA)-YNn  I  +  l)  )/CGSF(0CEL)  ) 

717  ODI(IDD)  =  CC+  IXNI(I  +  1)-  XNI  (lAA)  -(YNI(IAA)  -  YNI(Ul))* 

1  TANF(DDEL)  )*SINF  (BBE  ) 

WAXIE  =  DDI ( 1) 

,  ^  00  718  100  =  2,10 

718  WAXIE  =  MAXIF  (WAXIE,  001(100)) 

,ISP^r2»502,  BEJ,SGK,B2L,EPM,AKN, STAG, SOLID, WAXIE 
302  FORMAT  ( 1  Hi //52X 15HC ASCAOE  RESULTS// 

1  18H  DESIGN  PARAMETERS/ 

2  9H  ^BETA  *F6.2,  5H  DEG.  6X7HSIGMA  ■ 

3  F6.2,  9X  9HEPSIL0N  =  F5. 3,9X7HKAPPA 
u  17h  other  parameters/ 

5  IIH  STAGGER«F6.2,5H  DEG.  4X  lOHSOLIDITY  «F6.3  ,7X17HMAX.  THI 
6CKNESS  =  F5.3) 

IF  (NDPR)  140,140,141 

140  WRITE  OUTPUT  TAPE  2 , 304, XB, YB, PHEB, XT, YT, AHTO, TETD, PHT , XL , YL, AHLO, 
1  TELD,PHL,XP,YP,XQ,YQ,TAUO,XZ,YZ,AP.BP.BEPD 

3040F0RMAT  (3X4HXB  =F8 .4 , 1 IX 4H YB  *F8.4,  11X8HPHI(8)  «F8.4,/3X 

1  4HXT  *  F8.4.11X4HYT  »F8.4,  1 1 Xl 0HALPHA( T )  «  F6.2.  5H  DEG. 

2  2X10HTHETA (T)  =  F6.2,  5H  DEG.  2X8HPHI(T)  *  F8.4,/3X 

3  4HXL  =  F8.4.  11X4HYL  »  F8.4,  1  1X10HALPHA( L)  *  F6.2,  5H  OEG. 

4  2X10HTHETA(L)  =  F6.2,  5H  DEG.  2X  BHPHKL)  «  F8.4,/3X 

5  4HXP  =F8.4. 1 1X4HYP  =F 8 . 4 , 1 1 X 4HXQ  «F8 .4, T 1 X4HYQ  »F8.4,llX 

6  5HTAU  =F13.2,  5H  0EG.,/3X 

7  4HXZ  =  F8.4,  11X4HYZ  »  F8.4,  11X4HAP  «  F8.4,  11X4HBP  *  F8.4,11X 

8  12HBETA  PRIME  =  F6.2,  5H  OEG.) 

141  WRITE  OUTPUT  TAPE  2,303 

3030F0RMAT  (//42X  13HAIRF0IL  SHAPE  32X21HVEL0CITY  DISTRIBUTION/ 

1  88X  19HAT  DESIGN  CONDITION/ 

2  .10X1HI17X1HX19X1HY14X1  IHSLOPE,  DEG. 1 1 X6H( V/ V2 )  1 2X9H( V/ V2 ) **2//) 
1  =  1 

DO  90  1*1,42 

PRX=(XI I( I )-XI 1(41 )) /CHORD 
PRY=(ETI(  I)-ETJ(41))  /CHORD 
XNI  (  I  )  =  PRX*CNA+PRY*SNA 
YNI (  I  )  =  PRY*CNA-PRX*SNA 
DELI (I)=OELI (I )-AROT 
VVI ( I )=CV*VVI( I ) 

V2I (  I)  =  VVI ( I )*VVI ( I) 

90  CONTINUE 

WRITE0UTPUTTAPE2,305, ( I,XNI ( I ) , YNI ( I ) , DEL  I ( I ) , VVI ( I ) ,V2 1 ( I ) , 
11=1,42) 

305  FORMAT  (Ill,  2F20.4,  F20.2,  2F20.4) 

PARTIV  PARAMETER  CHANGES 
48  AKN=AKN+AKI 
CONTINUE 
EPM=EPM+EBI 
CONTINUE 
BZL=BZL+BZI 
CONTINUE 
SGK=SGK+SGI 
CONTINUE 
8EJ=BEJ+BEI 
CONTINUE 
ID=ID+1 

IF(IB-ID)  202,203,202 

PART^^V  ITERATION  SUBROUTINE  FOR  POINTS  ON  NEAR  CIRCLE 

SUBROUTINE  SSR  ( X, Y, VX, V Y, V ,C22 ,C23, A, B, PHN, PHEB, ITS , SNAL B,CNALB , 
ISIOE, ITSIDE,A2,  NEXPR) 

ITS  =  0 
ITSIDE=0 

IF(NEXPR)  8,8,517 

F0RMAT?//5H^  ITS^4X^HX  9X1HY  9X2HR1  7X3HAJ1  8X2HR2  7X3HAJ28X1HA 
1  9X1HB  8X2HVX  8X2HVY  8X2HA2  5X4HSI0E/) 


41 

40 

29 


203 


517 

510 


i 


80 


SOS 

509 

507 

9 

10 


8  C2=SINH(X) 

C3*C0SH( X ) 

C6=S INF(Y) 

C7=C0SF(Y) 

A2=C3-C7 
R1*C2/(A2*2. ) 

AJ1=C6/(A2«{“2. ) ) 

PX=C3*C7-1. 

PY*C2*C6 
PZ=PX*PX+PY*PY 
R2=PX/(-2.*PZ ) 

AJ2=PY/(2.*PZ) 

PH  =  X*C23+Y»C22+A*R1-B*AJ  1 
PS»-X*C22+Y*C23  +  A»AJ1  +  B»  R1 
OPH*PH-PHN 
VX=C23+A*R2-E*AJ2 
VY=C22-A*AJ2-B*R2 
V=-SIDE*SQRTF( VX»VX+VY*VV) 

A2*( ABSF(OPH)+AESF(PS) )/PHEB 
ITS=ITS-*-l 

IF(NEXPR)  507,507,508 

2,509,  ITS,X.Y,R1,AJ1,R2,AJ2,A,B,VX,VY,A2,  ITSIDE 
FORMAT  (15,  lOFlO.M,  E10.4,  15) 

IF(ITS-25)  9,100,100 
IF  I.000001-A2)  10,100,100 
DX=( VX*OPH-VY*PS) /(V*V) 

0Y=( VY*OPH+VX*PS)/(V*V) 

X=X-DX 
Y=Y-OY 
GO  TO  8 

100  YY=Y*CNALB-X«SNALB 

IF  (YY*SIOE) 101,11,11 

101  I TSIOE=ITS IOE+1 
IF(ITSI0E-3)  102,11,11 

102  YY=-YY 
XX=X*CNALB+Y*SNAL8 
X=XX*CNALB-YY*SNALB 
Y=YY*CNALB+XX*SNALB 
GO  TO  8 

11  RETURN 
END 

PART  VI  ITERATION  SUBROUTINE  FOR  INITIAL  STAGNATION  POINTS  ON  NEAR  CIRCLl 
SUBROUTINE  BAKER  (  SGK  ,XB , YB ,RB, ALB , A ,B ,R 1 , A J 1 ,R2 , A J2, ERH , A2 , 

1  I T,C22,C23,NEXPR) 

IF{NEXPR)515,515,500 

500  WRITE  OUTPUT  TAPE  2,501  ,  ,  ^  , 

501  FORMAT (//5H  IT  4X2HXB  8X2HYB  8X2HR1  7X3HAJ 1  8X2HR2  7X3HAJ2  8X1HA 
1  9X1HB  9X2HPS  8X2H0P  8X2HA2/) 

515  'CONTINUE  . 

IT=0 

PH£B=SGK*3. 141593 
SNALB=  SINF(ALB) 

CNALB=  COSF(ALB) 

XB=RB*CNALB 
YB=RB*SNALB 
4  C2=SINH(XB) 

C3=C0SH(XB) 

C6=SINF(YB) 

C7*C0SF(YB) 

A2=C3-C7 
R1=C2  /(A2*2.) 

AJ1=C6/(A2* (-2. ) ) 

PX=C3*C7-1 . 

PY=C2*C6 
PZ=PX*PX+PY*PY 
R2=PX/( PZ* (-2.  ) ) 

AJ2=PY/(PZ*2.) 

QX=A J2*C22-R2*C23 
QY=AJ2*C23+R2*C22 
(3Z*R2*R2+AJ2*AJ2 


81 


A=QX/02 

B=QY/CZ 

PH=XE*C23+YB*C22+A*R 1-B* AJl 
PS=-XB*C22+YB*C23  +  A*AJUB*R1 
OP=PH-PHeB 

A2={ ABSF{DP)+ABSF(PS) )/PHEB 
IT=IT+1 

IFCNEXPR)  51*4,514,503 

503  WRITE  OUTPUT  TAPE  2 , 504, I T, X8, YB ,R1 , A J1 ,R2, A J2 , A ,B ,PS ,0P,  A2 

504  FORMAT  { 1 5 , 1 CF 1 0 . 4 , 1 E 10. 4 ) 

514  CONTINUE 

IF{IT-25)  5.7,7 

5  IF{ERM-A2)  6,7,7 

6  A2»C3-C7 

R1X=(  1  .-C3*C7)/(A2*A2*2.  ) 

R  1Y={C2*C6)/ (A2*A2*(-2.)  ) 

AJ 1X=-R1Y 

AJl Y=R1 X 

PXX=C2*C7 

PXY=-C3*C6 

PYX=-PXY 

PYY=PXX 

PZX=2.*( PX*PXX+PY*PYX) 

PZY=2.* (PX*PXY+PY*PYY) 

R2X= (PZ*PXX-PX*PZX)/( (PZ*PZ)*(-2.) ) 

R2Y= {PZ*PXY-PX*PZY)/(PZ*P2*(-2. ) ) 

AJ2X=-R2Y 

AJ2Y=R2X 

QXX=-R2X*C23-4AJ2X*C22 

QXY=-R2Y*C23+AJ2Y*C22 

QYX=QXY 

QYY=-QXX 

0ZX=2.* (R2»R2X+AJ2*AJ2X) 

QZY=2.*{R2«R2Y+AJ2*AJ2Y) 

AX=( QZ*QXX-QX*QZX) /{ QZ*QZ) 

AY=(QZ*QXY-QX*QZY)/{ CZ*QZ) 

BX= ( OZ*QYX-QY*QZX)/(QZ»QZ) 

BY= ( QZ»QYY-QY*QZY) /(QZ*QZ 1 
PHX=C23+A*R1 X+AX*R1-B*AJ 1X-BX*AJ1 
PHY=C22+A*R1Y+AY»R1-B*AJ1Y-BY»AJ1- 
PSX=-C22+A*AJ1X+AX*AJ1+B*R1X+BX*R1 
PSY=C23+A*AJ 1Y+AY*AJ1+B*R1Y+BY*R1 
A2“PHX*PSY-PSX*PHY 
DX={PSY*DP  -PHY*PS)/A2 
DY= ( PHX*PS-PSX*OP  )/A2 
XB=XB-DX 
YB=Y8-DY 
^GO  TO  4 

7  RE=SQRTF{XB*XB+YB*YB) 

ALB=ATANF( YB/XB) 

IF(XB)28,45,45 

28  ALB=ALB+3. 141593 
45  CONTINUE 
RETURN 
END 
END 


5  6  1  5  45.  2.5  .5  .001  .01 


.4 


15  0  .245  1  0 


.  .END 
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32  J«1|NBE 

:=  (3. 141593/180. J  *  BEJ 


part  I A  CONSTANTS  FOR  NEAR  CIRCLE 
BE J*BE1 
002 

BE=  _ _ 

SNBE=SINF(BE) 

CNBE>COSF(BE) 

SGK=SG1 
003  K«1»NSG 

740  WRITE  OUTPUT  TAPE  2,  741.  BEJ. 

741  FORMAT  ( 1H1/26HCASCA0E  rIsULT^ 
1F6.3.//75H  B  ZERO  KAPPA 
2SS  BL  parameter/) 


SGK 

FOR  BETA  *F6.2.5H  DEG. 4X7HS IGMA  « 
EPSILON  SOLIDITY  THICKNE 


718 

751 

760 

755 


00  718  IDO  =  2. 10 

WAXIE  »  MAX1F  (WAXIEt  001(100)) 

PRAM=  V2I (1) 

00  760  1=  1.  19 
PRAM  »  MAXIF  (PRAMi  V2I(I)) 
WRITE 
FORMAT 


n«A  I  r  iri\F4nf 

OUTPUT  TAPE  2.  755.  BZL»  AKN, 
(3X,F6.2.2F10.3,2F14.3.F16.4) 


EPM«  SOLID*  WAXIE.  PRAM 


C  PARTIV  PARAMETER  CHANGES 


MW 


87  CONTINUE 

TNBE=SNBE/CNBE 

TNBE2*TNBE-BZL/2. 

CV2*(  1.+TNBE*TNBE )/I l.*TNBE2*TNBE2) 
CV*SQRTF(CV2) 

.  AKN*AK1 
00  U1  N>1«NAK 
WAXIT  =  0 
THIK  *  .04 
EPM-EP1 

IF(EPM)  99,98f99 

98  IF(N-l)  99t99tU1 

99  CONTINUE _ _ 


1  TANF(DOEL)  )*SINF  ( BBE  I 
HAXIE  =  DDI n ) 

DO  71  8  IDO  =  2t  10 
WAXIE  =  HAXIF  IWAXIEt  DDIIIDD)) 

IF  (WAXSG)  757,  757.  751 
IF  (THIK  -  WAXIE)  761,  751,  763 
IF  (WAXIT)  764,  764,  765 
EPM  =  EPM  -  .01 
GO  TO  99 

.01*1  (WAXIE  -  THIK)/(WAXIE  -  WAXIT)) 
WAXSG  =  100.0 
GO  TO  99 
EPM  *  EPM  +  .01 
WAXIT  =  WAXIE 
GO  TO  99 
PRAM=  V2I ( 1 ) 

DO  760  1=  1,  19 
PRAM  =  MAXIF  (PRAM.  V2I(I)) 

WRITE  OUTPUT  TAPE  2,  755,  BZL,  AKN, 

FORMAT  (3X,F6.2,2F10.3,2F14.3,F16.4) 

EPM=EPM+EBI 
THIK  =  THIK  ♦  .04 
WAXIT  =  0 
WAXSG  =  0 

IF  (THIK  -  .25)  99,  99,  48 
PARTIV  PARAMETER  CHANGES 
48  AKN=AKN+AKI 
41  CONTINUE 
^  BZL*BZL+BZI 
29  CONTINUE 
SGK=SGK+SGI 
3  CONTINUE 
BEJ=BEJ+BEI 
2  CONTINUE 
I0»ID+T 

^  ,  IF(IB-IO)  203,203,202 
.203_END  _ _ ! _ 


718 

756 

757 
761 

764 

765 


763 


751 

760 

755 


AKN,  EPM,  SOLID,  WAXIE, 


PRAM 
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